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Temperatures of Thermocouple Reference Junctions in an 


Ice Bath 


Frank R. 


(January 5, 


Errors obtained when using several sizes of ISA Type 
14 to 26 AWG 
Variables conside 
depth of immersion 
and type of insulation on the copper lead wires. 
Chromel 


AWG and of copper lead wires 
and maintained ice bath 
diameters and materials, 
june tion tubes 
the maximum copper and the 
normally used thermocouple wires. 

Depth of immersion 
degree of accuracy is to be 
type of insulation used a 


are given. 


minimum 

A small 
and diameter of the 
attained. 


1. Introduction 


A thermocouple used to measure temperature is, 


from its nature, a differential sensor which actually 
measures the thermal emf generated a result of 
the difference in temperature between its two 
junctions. It follows then, that the temperature of 
one of its junctions must be known if the measured 
emf is to be related to the temperature being 
measured. ‘This junction is known variously as the 
“cold,” “ice,” “reference” junction; and the fact 
that in measurements at low temperatures it may 
be the warmest part of interest in the measuring 
circult seems often to make no difference; this is 
still the “cold junction”’ to many, in spite of efforts 
of purists to have it called the “reference,” or when 
applic ‘able the “ice”? junction, The other is the 

‘measuring’ junction that often similarly is called 
the “hot” junction. In the discussions that follow, 
the “reference” “ice” junction refers to the 
electrical connection between the thermocouple and 
copper wires leading to the measuring instrument. 
This connection is made by mercury in the bottom 
of a Pyrex tube immersed in a bath of water-ice 
slush at, hopefully, 0 °C (32 °F). This is a con- 
ventional bath, descriptions of which can be found 
in the literature [1, 2].'. A discussion of the makeup 
of the ice bath including effects on its temperature 
of using tap or distilled water with tap or distilled 
water ice is given by McElroy in reference [3]. 


or 


or 


or 


Figures in brackets indicate the literature references at the end of this paper 


Caldwell 


1965) 


K thermocouple wires (14 to 20 
at different immersions in a properly prepared 
red include in addition to the wire 
2 to 9 in.), solid cross section of the Pyrex 
The wires studied represent 
thermal conductivities of anv of the 
amount of data on platinum is included. 
copper wire extremely important if : 


are 1 high 


The solid cross section of the Pyre x jt unction tubes and 
re of lesser importance. 


He warns that the reference junction tubes should 
be located properly in the slush bath if errors 
serious 2 to 3 °C” are to be avoided. 

The studies reported herein followed discussions 
in Committee AE Temperature Measurement 
Sensing, of the Society of Automotive Engineers, 
that led to the recommended practice of reference 
[2]. These emphasized the desirability of evaluating 
the errors resulting from varying geometries of 
tubes and wires in the ice bath. Conversely, such 
an evaluation would provide information on the 
sizes of tubes and wires and depths of immersion 
that would lead to certain tolerable limits . errors 

The purpose of this paper is, therefore, to point 
out and give quantitative information on errors that 
may arise from misuse of thermocouple reference 
junctions in ice baths, or on the limiting conditions 
that will cause an error not to exceed the maximum 
tolerable. Figure 1 is a sketch of an ice bath of the 
type used in this investigation. 

The temptation to cover the entire field of con- 
ventional thermocouple wires has been withstood, 
and as a result, the data obtained are confined almost 
exclusively to the ISA type K thermocouple wires 
(Chromel P and Alumel in this case) and copper 
leads. A few observations were made with plati- 
num wire and copper leads. The data presented 
thus are applicable to thermocouple elements of the 
lowest conductivity normally used (type K) and 
also of the highest (copper). 

The assumption is made in presenting this paper 
that the reader is familiar with the basic principles 
of thermoelectric thermometry to the extent, at 
least, that he can intelligently apply them to practice. 
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Figure 1. Ice bath and reference junctions. 


2. Practical Considerations 


If one end of a wire is immersed deeper and deeper 
into a bath of uniform temperature, a depth is reached 
at which greater immersion produces no further 
measurable or appreciable change in temperature of 


the immersed end. When the wire is in a glass tube 
as considered here and illustrated in figure 1, the 
same condition is reached, but at a somewhat greater 
depth. In addition to the depth of immersion, the 
thermal conductivity of the wires and material of 
the tubes, cross sections of the wires and tubes, 
depth of, and immersion of wires in the mercury in 
the tubes, type of insulation on the thermocouple 
wires or copper lead wires, and contact of the wires 
with the sides of the tubes all may affect the tempera- 
ture at the reference junctions. The condition of the 
ice bath also is important, and it should be renewed 
or replenished often enough that no measurable 
gradient exists in either the vertical or horizontal 
direction in the area of the tubes. 

Within limits, then, several choices are available 
for reducing the temperature of the reference junc- 
tion in an ice bath and consequently the error 
introduced by the difference between the temperature 
of the junction and 0 °C. The choices examined in 
this work are: cross sections of wires, of both copper 
and thermocouple materials; type of insulation on 
the copper wire leads; immersion of junctions in the 
ice bath; and cross section of the glass tubes. 


3. Materials and Procedure 


The first thermocouple wires examined were four 
sizes of Chromel and Alumel (ISA type K): 8-, 14-, 
20-, and 26-gage (American Wire Gage). Copper 
wires used with these were 14-, 20-, and 26-gage. 
The type K wires were oxidized, and the copper was 
enamel-insulated. Each increment between these 
wire sizes starting with the smallest represents an 
increase of about twice the diameter and four times 
the cross section. The dimensions of the several 
gages of wires used are given in table 1. 


TABLE 1. Diameters and cross-sectional areas at 20 °C of the 


wires used ® 


American wire 


gage 


Diameter Cross section 


m 
Zz 0. 1285 
14 06408 
20 . 03196 
26 01594 


mim 

3. 264 

1. 628 
S118 


. 4049 


mm? 
8. 366 
2. 081 
0. 5176 
. 1288 


in 
0.01297 

003225 

0008023 
. 0001996 


® From Smithsonian Physical Tables 


Smithsonian Institution, Washington, 
D.C., 1954) ninth revised edition; W. E 


Forsythe, Editor. 


Thermal conductivities of Chromel, Alumel, cop- 
per, and Pyrex No. 7740 are given in table 2. Here 
copper is seen to have over 21 times the thermal 
conductivity of Chromel P, about 13 times that of 
Alumel, 5.5 times that of platinum, and nearly 350 
times the thermal conductivity of Pyrex. 


TaBLE 2. Thermal conductivities at 0 °C of copper, Chromel 
P, Alumel, and Pyrex, and of copper and Pyrex relative to 
the others 


Material kpyrex/k 


Copper *(). 912 
Chromel P- 

Alumel 

Platinum 


Pyrex 7740 
Pyrex 7740 


0. 00286 
. 0621 
0384 
0157 


«C.8. Smith, The Physical Constants of Copper, Metals Handbook p. 1380 
(The American Society for Metals, 1939 

b Extrapolated from table in Hoskins Manufacturing Company Catalog 
M-61, C-A, 1961, p. 5, and converted from watts/em °C. Chromel and Alumel 
are registered trade-marks of the Hoskins Mfg. Co 

¢ Calculated from equation of Holm and Stérmer, Measurement of the thermal 
conductivity of a platinum specimen in the temperature range 19-1020 °C, Wiss. 
Veroffentlich. Siemens-Konzerns 9, part II, 312 (1930 

4 Private communication from T. W. Watson, Building Research Division, 
National Bureau of Standards, and converted from milliwatts/em °C. Pyrex 
is a registered trade-mark of the Corning Glass Works. This is the value that 
is compared with those of the alloy and copper wires 

e Smithsonian Physical Tables, ninth revised edition; 


Note 


table 555, p. 534. 


k is thermal conductivity in cal/em sec °C 


Thermal conductivities of all of the commonly 
used thermocouple materials are within the range 
between copper and Chromel. For wires of a 
particular diameter therefore, copper, having the 
highest conductivity, may be expected to cause the 
greatest error; and use of the low-conductivity 
Chromel alone, if possible, would result in the 
smallest. This latter condition is not practical, 
because copper is used almost universally for the 
lead from the ice bath regardless of the type of ther- 
mocouple used. For practical applications of copper 





in conjunction with other thermocouple wires whose 
combined thermal conduction lies between those of 
copper and Chromel, the errors also will lie between 
those resulting from the use of copper and Chromel 
alone. 

The diameters and solid cross sections of the Pyrex 
tubes used are given in table 3. Sizes of tubes are 
identified in the table and in the text by the letters 
of column 1. 


Diameters and solid cross-sectional the 


Pyrex tubes used 


TABLE 3. areas of 


0.D. L.D. 


Solid cross section 


mim 
22. 6 
17.1 
16.1 
11.4 


&.R4 
4.54 


Thermal electromotive force (emf) was measured 


with a type K—3 potentiometer and a galvanometer 


with a sensitivity of about 5 mm/xV. 

The ice bath illustrated in figure 1 was made in a 
large Dewar flask with a clear pli istic cover, drilled to 
receive the Pyrex reference junction tubes. This 


cover Was substituted for the often used cork stopper 


in order that the depth of immersion of the tubes 
could be measured more accurately. The intended 
depth of mercury in the tubes was ¥ in. when the 
wires were immersed in it. Mercury does not, how- 
ever, normally wet the Chromel and Alumel wires 
and so the depth of the mercury often was irregular 
across the tubes, deeper than in. on one side and 
not so deep on the other. This effect was greater with 
the larger wires, and undoubtedly in some cases 
affected the indicated temperature of the junction 
between the alloy and copper wires. Although 
only two tubes are shown in the cross-sectional view 
of figure 1, observations often were taken with 
several in the ice bath at once. 

All of the results presented herein were obtained in 
ice baths that were made and maintained according 
to the directions given in references [1] and [2]. 
Clear shaved ice and tap water were used, and 
water-ice slush filled the Dewar flask at all times 
during observations. This latter precaution is im- 
portant in the use of an ice bath, because the water 
at the bottom of a bath in which the ice is merely 
floating may be as high as 4 °C, the temperature of 
maximum density of water. This is not just aca- 
demic; it has been observed in many instances in 
which long reference junction tubes have been used 
in an inadequately maintained ice bath. The 
junctions in these long tubes that extended through 
the slush into the water below actually have been 
found to be at temperatures well above 0 °C, and 
replenishing the ice bath has corrected conditions. 

Use of tap water may, as discussed in reference [3], 
have some effect on the temperature of the ice bath, 


| 0.0137 
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and this should be considered in cases where extreme 
accuracy is required. For the purposes of this 
paper though, this effect is unimportant; because all 
observations are differential. The chief requirement 


here has been that the ice bath be at as nearly a 
uniform temperature as possible throughout. 


The procedure used was to immerse a tube con- 
taining one junction of, for example, a Chromel 
versus copper thermocouple so deeply in an ice bath 
that further immersion caused no detectible change 
in the thermal emf. This depth was found for the 
largest wires, both alloy and copper, and the largest 
glass tubes used, to be 9 in. or less: consequently, “the 
stationary leg of the thermocouple in all cases was 
immersed to a depth at which the top of the mercury 
in the reference junction tube was at least 9 in. below 
the top of the ice bath. The other leg was immersed 
initially to a depth of 2 in. above the mercury and, 
with the thermocouple connected to the potentiom- 
eter by copper leads, a reading of the thermal emf 
was taken. The two junctions and the circuitry are 
shown in figure 1. The immersion was increased by 
an inch and another reading was taken. This was 
repeated for immersions up to 9 in., or to depths at 
which the emf became zero or constant. In the 
latter case the constant emf, normally a fraction of a 
microvolt was taken as an indication of inhomo- 
geneity of the wires, and not a true reflection of a 
temperature difference between the two junctions. 
In such cases the constant readings were corrected to 
zero, and the same correction was added to the read- 
ings at all depths. 

The first observations taken were obtained by 
repeating the procedures just described for all 
combinations of gages of Chromel P and copper and 
of Alumel and copper. At least two sets of data 
were taken in this phase for each combination, and 
in one case nine sets were taken. Three sizes of 
Pyrex tubes were used for the variable-immersion 
leg; the smallest had a solid cross section of about 
71 percent of that of the largest. These are the 
first three sizes of table 3. 

Although the thermal conductivity of the Pyrex 
reference junction tubes is very low relative to that 
of Chromel or Alumel, and extremely low relative 
to copper, separate data were taken to determine 
the effect, if any, of a change in the solid cross 
section of the tubes on the temperature of the 
reference junction. 

Another set was taken to find whether or not the 
type of insulation on the copper lead wires to the 
reference junction had a measurable effect. In one 
case the single silk enameled insulation was used with 
22-cage Alumel, and in the other the silk was re- 
moved, leaving only the enamel. The same leads 
were used with a platinum wire 0.020 in. in diameter 
(24-gage), a commonly used size. The Pyrex tube 
in each case was E with a solid cross section of 
in. and 0.079 in. I1.D. This size was chosen 
as a compromise between the fragility of the smaller, 
thinner-walled tubes and the larger, more durable 
tubes of higher thermal conductance. 





4. Results and Discussion 


Data taken to find the effects of changing size 
ind immersion of the type K and copper wires are 
presented in figures 2 and 3. Plotted in each figure 
are curves of emf of one type K wire of one gage 
against copper of three gages versus depth of im- 
mersion of the reference junctions in the ice bath. 
Figure 2 is for Chromel versus copper, and the data 
for Alumel are given in figure 3. 

Chromel is thermoelectrically positive to copper, 
and Alumel and platinum are negative. This means, 
for example, that the positive Chromel element of a 
Chromel versus copper thermocouple must be con- 
nected to the plus terminal of a measuring instru- 


ment in order to get a positive reading when the 
measuring junction is at a temperature higher than 
that of the reference junction. Similarly, in the 
case of a copper versus Alumel thermocouple, the 
copper must be connected to the terminal of the 
measuring instrument. 

In the discussions and figures that follow all 
thermoelectric data are positive, i.e., the thermo- 
couples considered are Chromel versus copper, cop- 
per versus Alumel, and copper versus platinum. 
The thermoelectric power at 0 °C of Chromel to 
copper is 19.57 wV/°C (10.87 pV/°F), of copper to 
Alumel is 19.71 wV/°C (10.95 wV/°F), and of copper 
to platinum is 6.15 wV/°C (3.42 wV/°F). The fol- 
lowing list of abbreviations are those used in figures 
2,3, 4, and 5, and in the tables. 
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FiGuRE 2. Thermal emf of Chromel versus 14-, 20-, and 26-gage copper at various depths of immersion in an ice bath. 
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Although the data are somewhat erratic, especially 
at the lowest immersions, some facts are immedi: itely 
apparent. As expected, deeper immersion is seen 
to result in smaller error, or deviation of the tem- 
perature of the junction from 0 °C (32 °F). The 
error appears generally to be slightly smaller for 
Chromel than for Alumel, though this is masked in 
some cases by the lack of precision of the measure- 
ments; otherwise the patterns are quite similar. 

The most potent variable is seen from the figures 
to be the size of the copper leads from the ice june- 


tions. The relatively mild changes of solid cross 
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FIGURE 4. 
Alumel 


Copper was 22-gage SE; Alumel 


section of the Pyrex reference junction tubes had no 
discernible effect. The largest copper wire used, 
14-gage, was responsible for an error at an immersion 
of 2 in. of about % °C for 8-gage Chromel and 1 °C 
(1.8 °F) for 8-gage Alumel. At the same immersion 
and for the S-gage tvpe K wires, 20-gage copper was 
responsible for errors of 0.2 °C or less when used 
with Chromel and up to 0.4 °C with Alumel. A 
depth of 4 in. is sufficient to reduce the error with 
either Chromel or Alumel of 14 gage or smaller and 
copper of 20 gage or smaller to 0.05 °C or less. With 
two exceptions where the error was 0.02 °C, im- 
mersion of the 26-gage copper wire to a depth of 
3 in. gave an error of 0.01 °C or less. The errors 
resulting from use of copper leads of the latter two 
sizes at an immersion of 6 in. always are less than 
0.01 °C 

The data of figure 4 were taken to determine the 
effect if any of the solid cross section of the Pyrex 
tubes on the temperature of the reference junction. 
The thermocouple used here was copper versus 
Alumel, both 22 gage; and the 22-gage copper was 
single silk, enamel insulated. Each curve of the 
chart is for a specific immersion of the top of the 
mercury in the tube into the ice bath. 

Although the emfs developed here are generally 
lower at immersions up to 4 in. than are the ap- 
proximately corresponding values of figure 3, the 
curves of figure 4 do point out the relative in- 
sensitivity of the junction temperature to the size 
of the Pyrex tubes used. In using the procedure 
described earlier to correct for any residual emf, the 
emf at 9 in. immersion was made to be 0 uV. The 
emfs at immersions from 5 to 8 in. as a result, then 
turned out to be from 0.1 to 0.3 uV for tube cross 
sections down to about 0.018 in.? The values for 
the two largest tubes are bracketed above by those 
for 8 in. immersion and for 6 in. immersion below. 
Practically, this latter emf is of such a magnitude 
that the spread undoubtedly is within the limits of 
precision of the observations; thus the error over an 


020 024 
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Effect of solid cross section of Pyrex tubes on thermal emf of copper versus 
at various depths of immersion in an ice bath. 


was 22-gage. 


immersion of from 5 to 8 
constant at 0.2 uV or This corresponds to 
about 0.01 °C. Another set of data similar to this, 
but with 26-gage Alumel showed the spread to be 
from 0 to about 0.5 wV for immersions of 5 to 9 in., or 
about 0.02 °C. The wider spread here, as in some 
other instances, may have been caused by inhomo- 
geneities introduced by greater cold working of the 
smaller wire. 

Data presented in figure 5 are intended to point 
out any difference that might result from the use of 
different types of insulation. Curve 1 (solid circles) 
is for bare 22 Alumel and 22-gage single silk, 


in. can be considered 


sO. 


22-gage 
enamel insulated copper in tube E, and Curve 2 
(open circles) is for the same wires with the silk 
insulation carefully removed. The difference in 
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FIGURE 5. Effects of insulation of copper leads on thermal 
emf of copper versus Alumel and copper versus platinum at 
various de pths of immersion in an icc bath. 


@, denotes 22-gage SE Cu, versus 22-gage 
curve 3, 

» denotes 22-gaze EN Cu, versus 22-gage 
curve 3. 


Alu in curve 1 and versus 24-gage Pt in 
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100 





thermal emf at an immersion of 2 in. corresponds to 
about 0.06 °C. Similar observations were taken 
with copper versus platinum thermocouples. The 
difference in this case, in the reverse direction, is so 
small as probably to be within the limit of precision. 
The one curve (3) drawn through the observations 
thus can be taken as an indication of errors to be 
encountered when using 24-gage platinum and a 
22-gage copper wire lead. The error corresponds to 
about 0.07 °C at an immersion of 2 in., 0.03 °C at 3 
in., and vanishes at 4 in. 


5. Conclusions 


The foregoing data point out that the errors 
caused by insufficient immersion in an ice bath, even 
under the extreme conditions of large wires and 
glass tubes, can be small enough that they will not 
be the primary causes of inaccuracy in many cases 
of engineering measurements. In those cases where 
the errors are of an intolerable magnitude, use of 
smaller wires and greater depth of immersion of the 
junction below the surface of the ice bath can reduce 
the error practically to the vanishing point. The 


size of the Pyrex tube used to contain the reference 
junction can be varied, within reasonable limits, 
without any appreciable effect on the temperature 
of the reference junction. 


The type of insulation on the copper wire, single 
silk enamel or enamel alone were seen in the case of 
Alumel apparently to have a measurable effect on 
the temperature of the reference junction, but in a 
similar test with platinum, no such result was found. 
These results are, therefore, somewhat inconclusive 
although it does seem reasonable that an effect of 
the type of insulation used on moderate sized wires 
should be seen. 
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Superimposed Birefractory Plates 
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A study is made of the relations among the parameters involved in the passage of a 


beam of plane-polarized light through a pair of birefractory plates. 


For a pair of plates 


having known phase lags and a known angle between their principal axes, the angle of in- 


clination for 


necessary 


a plane-polarized beam to emerge plane-polarized is determined. 


In addition, the conditions for producing circularly polarized light are developed and pre- 


sented in a simple and convenient form. 


It is shown how the ‘‘combination”’ 
metric half-shades. 


1. Introduction 


Many optical instruments employ polarized light 
and associated birefractory plates as _ essential 
elements in their operation. For example a pair 
of quarter-wave plates is used in photoelastic stress 
analysis in order to observe a system of isochromatic 
lines in a stressed plate irrespective of the azimuth 
of the plate [1, 2].!. Exact quarter-wave compen- 


sators are especially needed for devices such as the 
one developed by Goranson and Adams [3] (based 
on the classical method first fully described by 
Friedel {4]) which was used initially to measure 


birefringence in stressed glass, and later for deter- 
mining the effect of pressure on the release of 
birefringence in some transparent materials [5]. 
Jerrard [6] has developed a similar instrument. 
The Friedel method has also been adapted in 
constructing an interference microscope [7]. 

The full advantage of the extraordinarily precise 
Friedel methods is realized only when using a com- 
pensator the optical retardation of which is almost 
exactly one-quarter wavelength. There are various 
procedures for constructing precise quarter-wave 
compensators. One of these is to utilize a suitable 
combination of two birefractory plates of arbitrary 
phase lag [8, 2]. A consideration of the properties 
of such a two plate combination brings out some 
relations which are of great interest and also provides 
a method for constructing such compensators. 


2. Combining Two Plates 


When a beam of plane-polarized light passes 
through a birefractory plate the emerging beam, as 
is well known, will be found in general to be elipti- 
cally polarized. In treating this problem and the 
problem of passing the beam through two super- 
imposed plates we shall follow the development of 


* Retired. 
1 Figures in brackets indicate the literature references on page 113. 
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It turns out that the “combination” quarter-wave 
plate shows a striking similarity in several respects, 


to the ordinary quarver-wave plate. 


plates can be used to construct several novel polari- 


Goranson and Adams [3] using their notation with 
some additions (see table 1). For the passage of 
plane-polarized light through a single plate, the 
four quantities, ¢’, 6{, ®; and i, which are defined in 
table 1 are related by the equations ” 


sin : sin 2i, sin #,, 


2 i, cos 2’ cos 26}, 


cot cot 2e’ sin 26, 


PD, 


tan 26;=tan 21, cos ;. 


2 Equations (1), (2), and (3) above have been obtained from eqs (6), (7), and 
h they 


8) in the paper by Goranson and Adams [3] using the substitutions whi 
have indicated. Equation (4) above has been derived from eq (5) of Goranson 
ind Adams following their procedures and making appropriate substitutions 


TABLE 1. Definitions of the symbols used in the text * 


Plates I, I]—are the first and second doubly re- 
fracting plate in the path of a beam 
of plane-polarized light. 

are the vibrational directions of the faster 

and slower wave respectively in plate I. 
These directions form a set of mutually 
perpendicular reference axes. 

are the vibrational directions of the faster 

and slower wave respectively in plate 
Il. These directions form a second 
set of mutually perpendicular reference 
axes. 

y—is the angle which denotes the orientation of plate 

I with respect to plate I] and is measured from 
a; to a; or from y; to a}. y is always positive 
and is limited to the first quadrant, 1e., 0°< 
¥<90°. 


ay and 71 


, , 
Qe and Y2 


« The definitions given in this table correspond to those given by Goranson and 
Adams [3] with some additions. It may be noted that the definitions given here 
pertain to the special case where the principal directions of plates I and II have 
been chosen as the reference axes, and as a result are simpler than those of Goran- 
son and Adams in many instances, 
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¥.—is the value of y when plates I and II are super- 
posed to produce a compound quarter-wave 
plate. 


’? 


are the ellipticities of the light emergent from 
plates I and II respectively. The elliptic- 
ity is defined as the angle whose tangent 
is +(b/a) where a and 6 denote the major 
and minor axes respectively. By con- 
vention, with the light coming toward one, 
€ is positive when the light vector along 
the emergent elliptical vibration moves 
counterclockwise and negative when clock- 
wise. —45°<e€<45°. 


is the angle between the vibrational direction of 
the entering plane-polarized light and aj. 
90° <i, <+90°. 

is the angle between the vibrational direction of 

the entering plane polarized light and a; or 
y —90°<i,<+90°. 

are the values of 

the emergent 

polarized after 


and IT. 


are the values of i, and i, respectively when 
the emergent beam of light is circularly 
polarized after passing through plates | 
and IJ. [It is understood in this case that 
plates I and II have been superimposed at 
the angle, y,, to make a compound quarter- 
wave plate. | 


and i, respectively when 
veam of light is linearly 
passing through plates I 


i, 


refers to light which has passed through plate I 
and is the angle between the major axis of 
the elliptical vibration and —90°<9i< 
+-90°, 

refers to light which has passed through plate I 
and is the angle between the major axis of the 
elliptical vibration and aj or y. —90° <6) 
<+90°. 


refers to light which has passed through plates I 
and II and is the angle between the major 
axis of the elliptical vibration and a, or 72. 
—90° <6,’ <+90°. 


is the phase lag of 


given birefractory plate. 
<< 360°. 


, Py—are the phase lags for plates I and II re- 
spectively. 


When a beam of elliptically polarized light such 
as the one emerging from plate I passes through a 
second birefractory plate, the emerging beam, 
well known, is also elliptically polarized. For an 
elliptically polarized beam of a given ellipticity and 
azimuth entering a birefractory plate of a given 
phase lag and azimuth the ellipticity and azimuth 
of the emerging beam are uniquely determined. 
Goranson and Adams have shown that the following 
relations hold: 


as is 


sin 2e’’=sin 2e’ cos $),;+ cos 2¢’ sin 263 sin ®, (5) 


| that ¢’’ is zero, and that i becomes 1,, 


, ’ Aen ae... = 
tan 26,’=tan 26, cos , meg, sin, (6) 


| where again the parameters are defined in table 1 


These equations in conjunction with tables 1 and 
2 and with proper regard to the significance of the 
+ and signs, are adequate to determine the cor- 
rect and unique values of ¢’’ and 6,’ as follows: 
Where an alternative in sign occurs (and where it 
occurs in other equations throughout the paper) the 
sign chosen depends on which of the two principal 
directions of plate IT is chosen as the reference axis. 
If y is measured from a; the upper sign is used, if 
y is measured from y; the lower sign is used. 

When either the upper or lower sign is used in (5), 
it will be found that the solution for ¢’’ is unique 
because, by definition, 45° <e’’ << +45°. When 
either the upper or lower sign is used in eq (6) it 
will be found that there are still two solutions for 
6’, but physical considerations show that only one 
solution is valid. The correct solution for 6,’ may 
be found by referring to table 2. 


TABLE 2. Jamits on the correct values of 20,’ which may be 
used to find a unique solution for eq (@) 
cot &,=cot 2’ sin 


where 0°<.@,<180° if e >0° and 180°<@,< 360° if «<0 


3. Restoration to Plane-Polarization 


First we determine what are the relations between 
the variables in order for the beam, originally plane- 
polarized, to pass through the two plates and to 
emerge plane-polarized. 

Here «’’=0. Hence by ( 


tan 2e’= Ftan %,, sin 263. 


Applying (3) we have 


tan ®, sin 26; tan %,, sin 263. (8) 
Since the angle between the corresponding princi- 
pal directions of the two plates is y, it follows that 


, , 


6,=6,+-y. (9) 


Substituting (9) and (4) in (8), we have (noting 
the inclination 


that will result in emerging plane-polarized light) 
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cot 2i —( ~tan® 2y) cos Oy (10) From (11) and (9a) it follows that 
“71 tan ®;, ~~") sin 2p 
tan 26;’=tan y/cos or —cot y/cos ©, 
After the correct minus or plus value has been 
determined for use in eq (10), there will be two values | and, from (11a) and (10a) when ¢;=4,,=@ 
of i,,, differing by 90°, which will satisfy the equa- 
tion. These two solutions are both valid and can | 26, -—2i,, or —2i, +180°. (13a) 
be associated generally with the principal axes of a ‘ 
single birefractory plate. For each of these solutions 
there will be separate values of 6/, 63, and 63’. These Ane. Pies 
two sets of values will be separated from each other For this simplified case (when # and y are given), 
by 90° just as the i, values are. The following the procedure for determining i,, and 6%’ is to cal- 
relations are useful: From (6) and (7) we find culate i, from (10a) and 6;’ from (13a). If desired, 
we may obtain 6; from (12a) and 6; from (9a). 
tan 26;’=tan 26; (cos ®,,-+tan ®;, sin ©;;). The values of the parameters of interest are pre- 
sented in two illustrative examples in table 3. In 
Multiplying and dividing by cos #,, we obtain: both cases $;=0,,—=70°; in one example y has 
: been set arbitrarily at 30°, and the second example 
tan 26;’=tan 263/cos ;. ( refers to the situation where ¥,.—41.20° which is 
the angle needed to construct a composite quarter- 
From (4) and (10) we obtain | wave plate. (It will be shown in section 4 how the 
value of y, is found.) Figure 1 shows the relation- 
tan : 1 eo he oe of the azimuths for transmitting plane-polarized 
cot 26 =( en 2y ) — oy’ (12) | light with two superimposed plates when 4, 
. si ig $,,=b=70° and y,=41.20°. 





We are now in a position, having two plates with 
#, and #,,, and placed at angle y to determine what 


inclination i,, of the entering plane-polarized “ms Paste 3. Values of the parameters of interest in producing 
plane-polarized light and circularly polarized light when it 


will cause the emerging beam to be plane-polarized; is given that &;= y= 0= 70° 


and also to calculate 6, , the azimuth of the emerging Two exam ples are shown: (1) for an arbitrary y=30°, and, (2) for the required 

beam. Formally stated the problem is: Given ®,, pe=41.2 The values are given fer the situation where the a’-direction 
‘ . . of plate Il is chosen as the reference axis 

$,,, and y to calculate i,, so that the beam will emerge 

plane-polarized; and also to calculate 6)’. Solution, 

calculate i,, from (10), 6; from (12), 6)-from (9) and 


6,’ from (11 


larized light 


3.1. Special Case for $,—4,, 0; =20.68,—60.32 y =34.33°, 


In many instances it would be convenient to cilities cia 
superimpose two identical strips or slabs. Here we 
would have $;=%,,, which we may call simply ®. 
For this case eq (12) becomes 


20.60°, —69.40° 
10.67°, —79.33 
51.87 38.13) 


cot 26; cot y or cot (90—y) 


Circularly polarized light ‘or 62 is undefined 


1 


from which Yircularly-polarized | 6;=3.80°, —86.20 
light ¥(=y-)= 


, a 41.20 @5= 45.00°, — 45.00 
26, Y, y+180° or 90°—y, 90°—y+ 180°. 2a) : a 


The parameters are obtained from the following equations (particularly appli- 
Similarly from (10) cable in this simple form since j=): 
Plane-polarized light Circularly polarized light 


cot 2i,, cot y cos ® or tan y cos , 


10a 


and, from (9) and (12a) 
26,=y, ¥+180° or 90°+y, 90°+y+180°. 
Also (incidentally) from (8) with 4, 


Py; 


sin 26;/sin 20;= +1. 





, 


1 


ww EMERGING PLANE-POLARIZED LIGHT 
re ae 
a4 


- Eg” = 34.33° 


> Yoen 
—— D <> i =-=- 34. : 
| ra si ~Pi ” 





ae \— 1sT AZIMUTH FOR ENTERING PLANE- 
POLARIZED LIGHT WHERE THE EMERGING 
LIGHT WILL BE PLANE-POLARIZED 


6.87° 








2ND AZIMUTH FOR ENTERING PLANE- 
POLARIZED LIGHT WHERE THE EMERGING 
LIGHT WILL BE PLANE~POLARIZED 


Figure 1. Illustration of the azimuthal relationships in the transmission of plane-polarized 
light by two supe rimposed doubly refracting plates when ®=70° and V=V, 41.20%. 


and [B] that there are two mutually perpendicular entering azimuths for which plane-polarized light 
merge plane-polarized and that the vibration direction will be rotated through the angle (42’’-ip,) which is the 
th, 4 


For any combination of birefractory plates, there- | out that the sense of this rotation will be either 
fore, there can be found two mutually perpendicular 
azimuths for incident plane-polarized light such 
that the emergent light is also plane-polarized. 
In this respect such a combination is like a single , 
TS sone plate. The effect of the combination In table 4 are presented values of the parameters 
differs (in general), however, from that of a single | for emerging plane-polarized light when the two 
plate in that the combination rotates the vibration | plates have different #;=%,;=%, but with the 


direction through the angle (6,’—i,,). It turns 


positive (counterclockwise) or negative (clockwise) 
depending on whether the plates are combined so 
that y is measured from the a; or y; direction. 


arbitrary parameter, y, set at 30° for all cases. 
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TABLE 4. 


Values of ip, and other parameters of interest in producing plane-polarized from two plates having different % 
but with the arbitrary parameter y set at 30° 


a 
for all cases 


The values are given for the situation where the a’-direction of plate II is chosen as the reference axis 


45 
30 
—39. 2 
140. 8 
— 30. 0 
150. 
33. 7 
— 146. < 
30. 
—150. —150. 


—150.0 —150. 


The parameters are obtained from the following equations (particularly 


) 2ip,’--- 
 ——— oe 
 ) eae 

4) 204’. _-- 


4. Simple Method for Producing Circularly 
Polarized Light 


A simple method for getting y¥, and i,, for 
#, and #,, begins as follows: Suppose we let 6’; 
equal 45 Then, from with ¢’’=45°, sin 
(2¢’+-%,,)=1, and therefore 


given 


(5) 


45 + ,,/2. (14) 
[f i, is chosen to give e’ given by eq (14), then 6; can 
be calculated and y=y, can be calculated by eq (9) 
to give #;=45°. This combination of i, and y, will 
(by the conditions leading up to eq (14)) vield 
e’’ =45°,1.e., emergent laevorotary circularly polarized 
light. 

By (14)—for 
polarized light 


this case of emerging circularly 


sin 2e’ -COS Py (15) 


and 
(16) 


tan 2’ cot Pyy. 


From (1) and (15) we have 


sin 2i - COS Pr sin 
and from (3) and (16) 


sin 26; -ecot Py,/tan ®,. (18) 


and from (9) and the above assumption that @3;=45° 


(9b) 


It is important 
regard to eqs (17) 


to note the following points in 
(18), and (9b) and the ‘‘combina- 
tion’’ plate considered here: (1) The equations admit 
solutions only if 90°<,+-@,<270°. (2) Careful 
examination shows that only one pair of solutions is 
valid in each equation, although there are additional 
apparent solutions. (3) The two solutions for y, are 
complementary. One solution for y, is obtained if 


applicable in 





90 105 135 150 
30 30 < é 30 
—90. 0 K “ 


90. 0 
—30. 0 
150. 0 

90. 0 
—YO0. 0 

30. 0 


—150.0 - 150. 


simple form since @;=1 


eq (10a) 
eq (12a) 
eq (13a) 
eq (va 


a, is chosen as the reference axis, the second solution 
is obtained if y; is chosen as the reference axis. (4) 
As noted in section 3.1 there will be a rotation of the 
principal directions of vibration for this ‘‘combina- 
tion’’ plate, the sense of this rotation depends upon 
the choice of y,. (5) The equations were developed 
to produce left-handed circularly polarized light, 
l.e., €'’=45°. However, for a given ‘‘combination”’ 
plate, light incident at azimuth, i=i,+90° will 


emerge as right-handed circularly polarized light, 
Be. €* 45°. 


4.1. Special Case for $;—4,, 


For two similar plates we put ®, , and (17) 


becomes 


coe 


sin ai... -cot ® 17a 


and (18) becomes 


sin 26;= + cot? ®. 


Now from (9) sin 26; therefore from (18a) 


cos 2y; 


cos 2y, + eot- ®, 19a 
Then, from (17a) and (19a) we may 


readily cal- 
culate the impinging 


azimuth and the requisite sep- 
aration, y,, to give emerging circularly polarized 
light. Incidentally, it should be noted that a solu- 
tion for y, is obtained when # ranges from 45 to 135° 
but not from 0 to 45° or from 135 to 180° 
As a supplement to the example shown in figure 
figure 2 shows the relationships of the azimuths 
for producing circularly polarized light from two 
supe rimposed plates when ¢;=6,;—==70°. Figure 
3 shows the rotation of the principal planes of vibra- 
tion in a “‘combination” quarter-wave plate when 
@=70°. As shown, the rotation may be either posi- 
tive or negative depending on whether the upper or 
ower sign is chosen in eq (19a) when solving for y,. 
Table 5 gives values of y,, i.,, and i,, for various 
values of @ when two similar plates are combined to 
make a quarter-wave plate. Here it is evident that 
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1ST AZIMUTH FOR ENTERING PLANE- 
POLARIZED LIGHT TO GIVE LAEVOROTARY 
CIRCULARLY POLARIZED LIGHT, 

i = 10.67°. 
wey 

i =i + ¥ = 51.87° 
“Ce “c, 


(FOR THE CASE WHERE Y is 
MEASURED FROM &') TO GIVE 
A COMBINATION QUARTER-WAVE PLATE WITH 
LEFT ROTATION (6," - i >0°) 

=. oe 2 “Pe 


y «41,0 ° 
c 





2ND AZIMUTH FOR ENTERING PLANE- 
POLARIZED LIGHT TO GIVE DEXTRO- 
ROTARY CIRCULARLY POLARIZED LIGHT. 
i = - 79.33°. 


= - 38.13° 


2ND AZIMUTH FOR ENTERING PLANE- 

POLARIZED LIGHT TO GIVE DEXTRO- 

ROTARY CIRCULARLY POLARIZED. 

i = - 10.67°. 

“<<. 

i #2 »T @ 98,39* 

“Ts , 

¥. = 48.80° (FOR THE CASE WHERE ¥_ IS 
MEASURED FROM y,') TO GIVE 

A COMBINATION QUARTER-WAVE PLATE WITH 

RIGHT ROTATION (6,"" - i 0°) 





a. 
rr 
~ 


Piel 
x — 7 i. 


FIGuRE 2. 


In case [A], ¥=WV-. is me 
1s Measured [ro 
is shown in eact 
rotary circularly p 
light (fig. 3 


the choice of y, is fixed by the phase lags of the 
plates. A striking feature of the relation of i,, to i, 
is brought out by the table, namely that (i : 
appears to be always +45°. As a matter of fact in 
a pair of plates adjusted to transform plane-polarized 
light into circularly polarized light the difference 
between the two i’s is exactly 45°. That is, it may 


£1 —Ip,) 


asured from a)’ and rotates the principal vibration directions to the left. 

ind rotates the principal vibration directions to the right by 
case that there exist mutually perpendicular azimuths of incidence which yield laevo- and 
irized light respectively, which are at 45° to the azimuths for the production of plane-polarized 


“Ps * 


1ST AZIMUTH FOR ENTERING PLANE- 
POLARIZED LIGHT TO GIVE LAEVOROTARY 
CIRCULARLY POLARIZED LIGHT. 

= 79.33° 


i =i +¥ = - 51.87° 
=C2 a os 


Illustration of the azimuthal relationships in the production of circularly polarized 
light by two superimposed doubly refracting plates when ® 


iv”. 


In case [B], ¥=¥-=48.80 
see fig. 3 It 
dextro- 


in equal amount 


be readily shown that 
(20) 


cos? 21, 


_— 
sin? 21,, D1 


The proof of relation (20) may be found in 
appendix 1. In this respect the “combination” 
quarter-wave plate behaves like an ordinary single 
plate with 90° retardation. 
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FIGURE 3. Illustration of the 


I'wo sit 
In [B], ¥=¥ 
is the same 
Note that the 
light (fig. 2 


48.80) 


TABLE 5 Values of We, ic, i 
ms J. Wey bey 


I'he values are given for the 


oo 
0.840 


0 
180. 0 


The parameters are obtained from the following equations 


rotation 
quarter-wave plates when ® 


uations corresponding to the two solutions of eq 
und is measured from 2’. 


situation where the a’ 


70 
0. 3640 
1326 
34 
. 66 
2. 40 
.20 
142 
342 
—. 3918 
—68. 66 
111. 34 


particularly applicable in tl 


1) We-- 


of the 


The m 


1ST AZIMU 
POLARIZED 


LIGHT WILL 


2ND AZIMUTH 
LIGHT 


ED 


POLS 


LIGHT 


LIGHT 


2ND 


ED 


70 


19A) are shown: In [Al], v 


Pry re 8) 


100 
—(. 1762 
0312 
—10.2 
169.8 
RS. 4 
44.2 
1, 028 
—0. 1736 


. eq (19a 


. eq (17a 


eq (10a 


TH FOR 
LIGHT WHERE 
BE PLANE-POLARIZE 


WILL 


WILL 


IMU TH 


direction of Plate II is chosen 


lis simple form since 4, 


PLANE -POLARIZED 


ENTERING 


THE 


FOR ENTERING 


BE PLAN‘. -POLARI 


FOR 


LIGHT 


41.20 


wgnitude of the rotation is the same in both cases 
is that in figure 1 where the azimuthal relations of the principal vibration directions of the plates are shown 
izimuths in each case are mutually perpendicular and at 45 


is the reference 


110 
—(). 364 
. 1326 
—21.34 
158. 66 
82. 40 
41.20 
1. 142 
—(). 342 —0. 500 
3918 . 706 
66 54.8 
34 125.2 


144.8 
70. 6 
35.3 

1. 412 


$11): 


AXIS 


LIGH< 


principal vibration directions in combination 


and is measured from a2 


Case [A] 


to those for producing circularly polarized 


p,, and related parameters in producing circularly polarized light from two plates having different 
2; 


0 
180. 0 





In another respect also the “combination’’ plate 
resembles a simple quarter-wave plate, namely that 
when a beam of plane-polarized light is passed 
through, the azimuth of the ellipse of the emerging 
beam is independent of the inclination i, of the 
entering beam. For proof of this interesting relation 
see appendix 2. 

There is a third way in which the “combination” 
plate resembles an ordinary quarter-wave plate. 
To show this, we first recall that by (1) a beam of 
plane-polarized light passing through a plate with 
® equal to 90° emerges with ¢ equal to i. For our 
“combination” plate—let us put 

A=i,- (21) 
It turns out, as shown in appendix 3, that whereas 
with the ordinary quarter-wave plate e’=i, for the 
combination quarter-wave plate 
i, i,, A. 


(22) 


Thus, there is the same relation, between ¢’’ and A 
as there is between ¢’ and ifor the simple quarter-wave 
plate. 


5. Design of Half-Shade Systems 


Precision measurements in polarimetry usually 
depend upon making a brightness match between 
two parts of a field of view. Half-shade devices 
have been constructed which permit the observer 
to ascertain very precisely the azimuth of plane- 
polarized light, and the azimuth and ellipticity of 
elliptically polarized light [6, 9, 10]. As examples, 
the lippich half-nicol and the Cornu-Jellet split 
nicol are widely known as azimuth half-shades, and 
the Brace half-shade has been used extensively in 
the measurement of ellipticity. The superimposition 
of doubly refracting plates makes it possible to 
construct several novel half-shade devices. These 
half-shades, partic “ul: urly ads apted to monochromatic 
analysis, are interesting in themselves, and they 
serve to illustrate some interesting features of the 
“combination”’ plates discussed thus far. 


5.1. Azimuth Half-Shades 


The inherent rotation of the principal vibration 
directions of a ‘combination’? quarter-wave plate 
which was discussed in sections 3.1 and 4.1 may be 
used to construct an azimuth half-shade. It is 
convenient to denote by a; and yz the principal 
directions (after rotation of als aevorotary compound 
pli ite, and by ap and 9 YD the prine ipal direc ‘tions (after 
rotation of a dextrorots iry compound plate. 

(a) In figure 4 it is shown that two compound 
plates may be cut along yz and ys and these two 
directions placed in adjacent positions alone the 
vibration direction, ?, of a nicol prism. The com- 
bination of nicol and two plates in fixed position then 
acts effectively as an azimuth half-shade with half- 
shadow angle of 2.6;'’—i The characteristics of 


110 


this half-shade are analogous to those of the Lippich 
half-nicol or the Cornu-Jellet split nicol. As an 
equivalent arrangement to the one shown in figure 4, 
a, may be interchanged for YL; when, 
time, aj is interchanged for 75. 

(b) The half-shade described in section (a) above 
depends upon the rotation of the optic axes which is a 
unique feature of “Combination” plates, and the two 
parts of the biplate may have an arbitrary path 
difference. On the other hand, half-wave plates 
have long been known for their ability to rotate the 
plane of polarization of incident light [11]. The 
Laurent plate is a half-wave plate placed partly over 
a polarizing prism and rotated in fixed combination 
with the prism [11]. Chauvin has described an 
analyzer consisting of the same combination [9]. 
The advantage of these half-shade devices lies in the 
fact that the sensitivity may be varied by varying the 
angle between the vibration direction of the polariz- 
ing prism and one of the principal directions of the 
half-wave plate. 

A half-wave plate may be simply constructed by 
superposing two of the ‘“combination’’ quarter-wave 
plates described in this paper. However one of the 
quarter-wave plates must be dextrorotary and the 
other laevorotary and they must have equal amounts 
of rotation of the optic axes. The two plates must 
be combined so as to compensate for the inherent 
rotation of each one. 

(c) It is known that when a doubly refracting 
material is placed with its principal vibration direc- 
tion in a diagonal position between two crossed 
quarter-wave plates, that the combination acts as 
as an optically rotary material [3]. Two specimens 
may be cut from the same doubly refracting sheet 
and positioned so that the fast direction of one 
specimen adjoins the slow direction of the other. 


at the same 





o* 7 
a 


FicgurE 4. An azimuth half-shade consisting of two compound 
plates in fixed combination with an analyzing nicol prism. 


Incident plane-polarized light with azi 
the field of view to 
azimuth shown 
extinction, 


shown at D causes the top half of 
show tot 3] extinetio " Inei dent plane-polarized light with 
it L causes the bottom half of the field of view to show total 











TWO DOUBLY REFRACTING 
SINGLE PLATES 





JARTER-WAVE PLATE 





Fiat RE 5. An azimuth half-shade consisting of a divided 
doubly refracting material between two crossed quarter-wave 
plates. 


The analyzing 
two halves « 


nicol rotates inde pe ndently and is shown in a position where the 
f the field of view have equal brightness. 


If these two plates are placed so that their principal 
vibration directions are in a diagonal position be- 
tween two crossed quarter-wave plates, the field of 
view being evenly divided, this fixed combination 
may be employed as an azimuth half-shade whose 
half-shadow angle is equal to ®, the phase lag of the 
divided plate. One possible arrangement is shown 
in figure 5. 

This half-shade acts then in a manner analogous 
to two plates of quartz cut perpendicularly to the 
optic axis, one plate being laevorotary and the other 
being dextrorotary. The combination described 
above acts effectively as a Soliel or Nakamura 
half-shade [6]. 

The arrangement shown in figure 5 is correct 
when using quarter-wave plates with no unusual 
features. However, in constructing the half-shade 
from quarter-wave plates described in this paper, 
the inherent rotations of the plates themselves must 
be recognized and accounted for. For example, in 
the arrangement shown in figure 5, it is possible to 
use a compound laevorotary plate for the first 
quarter-wave plate, and to set the axes so that 
ay, is fixed at the azimuth of a’ shown in the figure. 
In this case a dextrorotary compound plate must 
be used for the second quarter-wave plate, y’ (not 
yp) of this second plate will still have the azimuth 
of y’ shown in the figure. 


5.2. An Elliptic Half-Shade 


The Bravais biplate is well known as an elliptic 
half-shade [10]. It consists of two quarter-wave 
plates with their edges in juxtaposition and the 
fast ray of one plate parallel to the slow ray of the 
other. This biplate serves to show when elliptically 
polarized light has been completely restored to 
plane-polarized light, because the two halves of the 
field match at all azimuths, the brightness level 
































Fiat RE 6. Bravais biplates. 


A An elliptic half-shade in which the two parts are single 


opposite order quarter-wave 
B_ An elliptic half-shade in which the two parts are compoun 
and opposite order (quarter-wave), 


alone varying, but the match is destroyed as soon 
as the incident light becomes elliptically polarized. 
A Bravais biplate is shown in figure 6A. 

Because the azimuth of this biplate may be varied 
with respect to the analyzing nicol, it is perfectly 
all right to use compound quarter-wave plates with 
their associated rotation of the principal axes. It is 
important, however, that both halves of the biplate 
have the same amount and sign of rotation. Figure 
6B shows a Bravais biplate constructed from two 
laevorotary quarter-wave plates. The biplate might 
as easily be constructed from two dextrorotary 
quarter-wave plates. 


5.3. Half-Shades Uniquely Adapted to the Friedel 
Method 


The third Friedel method has been employed by 
Goranson and Adams for the precise determination 


of optical path difference [3]. According to this 
method, plane-polarized light is made to pass through 
a doubly refracting specimen with one of the prin- 
cipal vibration directions of the specimen set at 45 
deg with 1 respect. to the vibration direction of the 
polarizer. The light is then made to pass through 
a quarter-wave plate with one of its principal direc- 
tions parallel to the vibration direction of the 
original plane-polarized light. The light emergent 
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FIGURE 7. 


A 
B 


There is no divided 
Two compound 


field. 


Friedel method III. 


The analyzing nicol must be rotated to the left for complete extinction 
quarter-wave plates produce a divided field 


The analyzing nicol must be rotated to the left to make 


1 brightness match between the two halves of the field. 


( A divided field is 


produced by an auxiliary doubly refracting plate. The analyzing nico] must be rotated to the left 


to make a brightness match between the two halves of the field 


from the quarter-wave plate is plane-polarized, and 
an extinction position may be found by rotating the 


analyzer. The measured angle between the azi- 
muths of the incident and emergent plane-polarized 
beams of light is equal to one-half the phase angle 
of the specimen. The arrangement of the optical 
elements for this method is shown in figure 7A.* 


Although no half-shade device is shown in figure 7A the original instrument 
designed by Goranson and Adams [3] was equipped with a Wright biquartz 
wedge plate [12]. 


One may employ a compound quarter-wave plate 
described in this paper instead of the more familiar 
type of quarter-wave plate which is usually obtained 
by splitting sheets of mica to the proper thickness. 
In this case, however, it is necessary to know the 
sign and amount of the inherent rotation of the 
quarter-wave plate itself and to account for it in 
setting the zero position of the analyzer. 

(a) An effective half-shade may be obtained if 
two compound quarter-wave plates (one dextro- 
rotary, and one laevorotary) are placed with their 
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edges in juxtaposition so as to divide the field and 
with the a’ direction of each alined with the principal 
vibration direction of the polarizer. This arrange- 
ment is shown is figure 7B. In this case there is no 
correction to be made to the zero position of the 
analyzer, and, in measuring, the analyzer is rotated 
until the two halves of the divided field have the 
same intensity. 

(b) Another half-shade device may be obtained if 
a second doubly refracting plate is placed after the 
specimen so that it covers half the field and has its 
principal vibration directions alined with those of 
the specimen. The auxiliary plate may be selected 
so that the rotation incurred by its phase lag just 
balances the inherent rotation of a compound 
quarter-wave plate. Actually, @ of this auxiliary 
plate must have the magnitude 4| 6;’—i,, 
6,’ and i,, refer to the 
late. A typical arrangement is shown in figure 7C 
= the case of a dextrorotary compound quarter- 
wave plate. In making a measurement, as with the 
half-shade device discussed in part (a), the analyzer 
is rotated until the two halves of the field have the 
same intensity. 


“combination” 
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7. Appendix 1. Relation Between i. and 
i,, in an ‘‘Adjusted’’ Combination Plate 


Given: and the 


(When y is measured from a;) 


the four quantities %, y,, i,, and i, 
three equations. 
t-cot* & 


cos 2y, (19a) 


tany,=—tan 2i,, cos & (10a) 


sin 2i,,=cot ®. (17a) 


To prove: 
(20) 





, where | 
quarter-wave 


from (19a), 
1—cot? ® 
1+cot? @ 


1—cos 2y, 
1+cos 2y, 


sin 1 2We 
cos QW. 


sin? 6—cos? ®- 


-tan? y, 


dividing by cos? ® gives 


] a tan’ y, 


cot? ® cos? ® 


is 
14 tan’ y, 


cos’ ® 


cot? &= 


From (23) and (10a) we have, 


1 — 
cot" b= —==COs* 21, (24) 


1+tan? 2i,, _ 


From (24) and (17a) we have sin? 2i, 
oe 


(The proof when y is measured from y; 


CC Ss? pi 


Pi* 


is similar.) 


8. Appendix 2 


Given eq (6) and other relations pertaining to the 
“combination”’ quarter-wave pl: ite, it is required to 
prove that the resulting azimuth, 6,’, is independent 
of the inclination, i,, of the entering plane-polarized 
beam. 

By (3) the eq (6) becomes 


ae tan®@sin#, . 
tan 26, =tan 20, cos 6— — sin ®. (6a) 
cos 26, 


In order 


follows: By 


to eliminate 6; 
(9) and (4), 


and 6%, 


we proceed as 
we have 


tan 2i, cos $+tan 2y, 


tan 26; . . 
1—tan 21, cos ® tan 2y, 


Also, 


we obtain from (9) 


sin 26; ] 


cos 26, cot 26, cos 2y,- 


-sin 2y, 
From this, with (4) and (19a) 
sin 26; 
cos 26; cot* & 


1—ec th 
tan 2 \ ais 


1, cos ® 


and, by (19a) 


tan 2y, (= = 


Now substitute tan 2y, from (28) in (25), and put 
the resulting value of tan 265 and also the value of 


\ 1—cot* ® 
cot? & 


(28) 


sin 2y, ) 
cos 2y, 
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sin 26)/cos 26; from (27) in (6a), and we obtain an 
expression equivalent to the second term of (6) or 
(6a). Multiply the numerator and denominator 
by sin? @ to get a modified expression, which we call 
A, Le., 
tan 26;’=A. 
It will be found that A cos @/,—cos 2@turns out 
to be simply unity (the quantities i, and ® happily 
dropping out Therefore, 


ae —cos 26 
tan 26, : tan 2y,. cos ®. (29) 


cos ® 


The azimuth of the elliptically polarized beam is 
therefore independent of i). 


9. Appendix 3. Ellipticity of a Light Beam 
Emerging From an ‘‘Adjusted’’ Combina- 
tion Plate 


We have $,=4,,=@ and y is measured from aj. 
By substituting (1) and (2) in (5) we obtain 
sin 2e’’ =sin 2i, sin ® cos ® 
-cos 2i, sin ® sin 265/cos 26}. (30) 
From (9), 


sin 26; sin 26; cos 2y+cos 26; sin 2y, 


(9c) 
and, introducing (9c) into (30), 
sin 2e’’=sin 2i, sin ® cos ® 


L_cos 2i, sin ® (tan 26; cos 2y+sin 2y). (31) 


From (19a), 


sin 2y,= y1—cot! %, (19b) 


and, introducing (4), (19a), and (19b) into (31), 


sin 2e’’=sin 2i; sin & cos ® 


+cos 2i, sin ® (tan 2i, cos ® cot? & 


y1—cot* ®). (32) 


noting that, in general, 


sin ® cos ® cot? b=cos' &/sin & 


sin ® cos ©-+-cos® &/sin P=cot & 


| 
| 


sin Pp \ 1—cot* @ \ l -cot? P, 


sin 2e’’ =sin 2i, cot + cos 2i,41—cot? &. (33 


From (24), 


+sin 21), y1—cot* ®. (24a) 
From the geometry the sign can be either plus or 


minus; taking the case where the sign is minus— 


. 6 tv 
sin 2e 


sin 2i, cos 2i,, 


cos 21, sin 2]),. 


(34) 


sin 2e’’ =sin (2i,—2i,,), (35) 
(22) 
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Accuracies obtainable in measuring resistance (conductance) at radio frequencies have 
lagged behind those obtainable in measuring reactance (susceptance), because frequency- 
dependent changes in the values of resistors are not as well known as are such changes in 
capacitors, and because no satisfactory method has existed for comparing resistors to capaci- 
tors at radio frequencies. An instrument based on the twin 7 null circuit can be made 
self-calibrating at a given frequency, allowing conductances in the proper range to be meas- 
ured directly in terms of a change in capacitance, without any need for accurate knowledge 
of the values of the circuit elements. Analysis of a practical circuit shows that such an 
instrument is potentially capable of measuring conductance to within a few parts in 10° at 
frequencies below about 15 MHz, and that accuracies actually obtained will be limited 
only by the accuracy to which changes in effective capacitance are known. Conductances 
between 0.01 and 0.0001 mho were measured at 107 radians/see with an estimated error of 
0.05 percent, which is the error associated with the measurement of changes in -effective 
capacitance at that frequency. 


1. Introduction 


[t has long been known that capacitors are superior to either inductors or resistors for use 
as impedance standards at high frequencies, because changes in their effective values with 
increasing frequency can be more accurately and more easily evaluated than can such changes 
in either resistors or inductors. For this reason, the High Frequency Impedance Group in 
the Electronic Calibration Center at NBS-Boulder, has decided to refer all its measurements of 
impedance at high frequencies to standards of capacitance [1].!_ While either null or resonance 
techniques can be used to compare capacitors and inductors to derived standards of capacitance, 
no satisfactory technique has existed for comparing resistors to capacitors at radio frequencies. 
The existing null techniques for measuring resistance require prior knowledge of the effective 
values of at least one resistor at the frequency at which the measurement is to be made. The 
resonance techniques, while theoretically capable of measuring resistance directly in terms of 
a change in capacitance, are at present capable of high accuracy only when measuring very 
high or very low values of resistance. Because of the lack of a satisfactory measurement 
technique, the accuracies obtainable in resistance measurements have lagged behind those 
obtainable in the measurement of reactance. 

An investigation of the twin 7 null circuit has shown that it is possible to construct an 


instrument capable of measuring resistance directly in terms of a change of capacitance, without 


prior knowledge of the values of any of the circuit elements. 

This paper describes an instrument which was constructed for this purpose, and presents 
results which indicate that resistances, in the proper range, can be measured to an accuracy 
which is limited only by the accuracy to which the change in capacitance is known. 


2. Circuit Equations 


2.1. Ideal Case (No Residual Impedances) 


The circuit to be considered (fig. 1) is adapted from one which was described by Tuttle 
|2], and applied to a commercial instrument by Sinclair [3], in 1940. The circuit is essentially 
the one which Woods [4] used in his VHF Dual Admittance Bridge in 1950. Tuttle has shown 


lards Engineering Division, NBS Boulder Laboratories, Boulder, Colo 
1 brackets indicate the literature references at the end of this paper 
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FIGURE 1. Twin T circuit using pure com- 
ponents 


TERMINAL 
PAIR B 


| 


= ae 
‘. i. | 


Se 


that the condition for zero transmission through the twin JT (null condition) is given by 


aa, tA 0. 


Ly sf T ; a 


l ; l ' : 
In the present case, Z;=—; Z,;—=~—. ete., and the real and fimaginary parts of (1) give the 
JOC; / 


two null equations: 


wR.CC3 (1 + 


As in the Woods Bridge, a measurement can be made by connecting an unknown admit- 
tance Y,—@,+jwC,, either to terminal-pair A, or to terminal-pair B. When Y, is connected 
to terminal-pair A, C, must be replaced by C)+C,, G. by @+G,, and C; by C%, in the null 
equations. Similar changes (indicated by double primes) must be made when J, is connected 
to terminal-pair B. Subtracting the initial null equations from the final ones yields the four 
expressions 


(4) 


where 





and subscripts A and B denote the terminal-pair to which Y, is connected. 
If P is defined as 

CC: 

P=R, 7 = 


4 


G.,=w0°PAC, 


—] ; 
G, ACh. 
B P 2 


If a conductance is measured first on terminal-pair A, then on terminal-pair B, 


G, 1 Gep 
and 
] 


w PAC3;=— > AC, 
5 P 2 


sat ‘ 1 
If P is adjustable so that it can be made equal to ~, 


—1 AC, 
w AC; 

and 
AC,=—AC; 


If P is adjusted to make the same change in capacitance necessary to null the instrument when 


G, is measured on terminal-pair A, as when it is measured on terminal-pair B, then P=-> and 


G,4=wAC; (10) 


ee 5 A (11) 
While P has not been eliminated from the calculation of G,, it is possible to evaluate it with 
great accuracy, allowing conductance to be measured with accuracies limited only by the 
accuracy to which a change in capacitance is known. 

It should be mentioned that a conductance may be measured in terms of a change in 
capacitance without making P equal to 1/w. If G,,=G@,,.=—G,, multiplying eq (8) by eq (9) 
vields the expression 


wAC,’AC;, 


and P is eliminated from the calculation of G,. Because the number of measurements needed 
to determine G, is doubled, and because accuracy decreases as the value of P departs from 
1/w, this method is less effective than the other, especially when more than a few measurements 
are to be made at a given frequency. 


2.2. Realistic Circuit (Including Residual Impedances) 


In general, residual impedances cause the effective values of the circuit elements to change 
with frequency and cause losses in the fixed capacitors, variable losses in the variable capaci- 
tors, and some reactance to be associated with the fixed resistors. 
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TERMINAL 
PAIR A 


FiGguRE 2. Twin T circuit using impure com- 
ponents. 


TERMINAL 
PAIR B 





The circuit of figure 2 is a more realistic representation of the actual circuit. In the 
figure, the value of each circuit element is taken as its effective value at the frequency at which 
the instrument is used. If FR, is lumped with G, for the present, and 2; with G;, the null 
equations for the circuit of figure 2 are 


l 


GG, ( G4 Gs 5 )— «(143+ Cy@,) (Cy+-Cs+-Cs) — eC; ( G,4+-G5+ : ) 


R, 


] 


| Gs) w (S4+0,4,)(C.— pt Ort-Cs )\—w?CC,(G,+ E244; 


R, 


W Ls 


and 


G,G;(C,+-C;+C,)4 y 4 C3G,)( G4 G7, 


C5 )+( nt! "Gy ) (G,+-G.4 


In the appendix, these equations are reduced to give the four expressions 
aw ( R,¢ 4 R,¢ ‘ (AC 9 + A AG), 


w*Ro(C.+C)(RiCi— ReC,) |— wR, C, + RsC,)(ACi+C.s)—AG;, 


( (RC, : RCs) | - [RiC, T RC RC, T RCs 1 R5(C; T C3) Gra, 


-AC é T [2z.0, 1 R Cc ie OP T RGF R2(C T C.)iG«. (17) 


For the instrument whose construction is described later, these equations are accurate to 
within a few parts in 10° If the instrument is to be self-calibrating for conductance 
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measurements, eqs (14) and (15) must be reduced to the form 

Gia t+AG,=w'PAC,, 
and 

G.2+AG;=— D 2, 
In order to accomplish this, define 

K=Rh,C,+ R,C,— RCA RC. 
Equations (16) and (17) then become 
ACJ1—w? R2(C2+ C2) (RyC,— R,C;) |—[K+R;(C5+ C3) Ges 19) 

and 


OC. —AC.+[K—Rx(C.+C;)]G,s- 20) 


When measuring a nearly pure conductance, sufficient accuracy is obtained for C,, if the term 
wPho(C,+C)) (RyC,— R,C,) is neglected. Then (19) becomes 


Cry —AC,—[|K+R;(C;+ C5) |Gea. (21) 


If a nearly pure conductance is measured on terminal-pair A, then on terminal-pair B, using 


the same capacitor operated over an identical range for both measurements, then 


G..=G.s=G 


z 
and 


Since Ces Cz, 
AC,—|K+ R(C+C’) |G, 


from which is obtained 


which can be evaluated. 
Now, R,C,, R3C3, and R,C, are likely to be at least an order of magnitude smaller than 
R,C;. When this is true, K=P,C,, 
(A+ R;(C3;+C) |G, 
and 
3 = [Ar R(C2 T C; IG zp. 


Equations (14) and (15) can now be rewritten 


G4 =o PAC;,— ov K[ K+ R;(C3+ CD |G.4—AGe 
and 


(1+ eK R(Cy+ C2) |— wo RC, + RoC) K— Re Cr+ C2) Gn AG. 


Grp 4 

If the terms wA|AK+&;(C;+C;)), &AR(C,+C)), and w(R,C,+ R;C;)[AK—R(C2+C)| are in- 
significant compared to unity, the expressions for the conductance reduce to the desired form. 
If AG,=AG;, or if both are insignificant compared to G,, P can be determined as it was in 

the ideal case. It is possible to insure that AC,—AC;, and that AG,—AG; in the two measure- 


119 





ments, by using the same variable capacitors operated over an identical range to measure G, 
and C, in both cases. This requires both G, and the capacitors which measure @, and C, to 
be switched from terminal to terminal. 

It is possible to adjust P to within the limits of stability and resettability of the capacitor 
used to measure G,. In a good two-terminal air capacitor equipped with a precision coaxial 
connector, the error is no more than a few parts in 10° of the maximum capacitance. Since 
changes in effective capacitance are not known to better than several parts in 10* at present, 
adjusting P in this way allows one to measure conductance to the same accuracy as that to 
which a change in effective capacitance is known. 

With the constants of the instrument known, it is only necessary to evaluate the residual 
impedances associated with the external variable capacitors, perhaps by the method of Field 
and Sinclair [5], in order to make quite good measurements of either conductance or reactance. 


3. Construction of the Instrument 


In converting the twin 7 circuit to a working instrument, several facts must be considered. 
Since it is necessary to move the standard capacitors from one pair of terminals to the other 
when evaluating the constants of the instrument, the standard capacitors cannot be made an 
integral part of the instrument. If the capacitors are not to be built in, more than one pair 
of terminals must be provided to allow paralleling the unknown admittance and the capacitor 
which measures C,. Multiple terminals are also desirable because they simplify calibration 
of the variable capacitors used with the instrument. If multiple terminals are to be used, it 
is desirable for each pair to connect to a common point, the “bridge corner,’ and for the 
conductors to be mutually perpendicular in order to reduce mutual impedances as much as 
possible. The fact that it is possible to connect external variable capacitors of different ranges 
for C, and C; requires changing the values of L and G, necessary to provide initial balance. 
Since multiple terminals are already necessary, it is just as well to provide for the attachment 
of external standards for Z and G,. As a result, the actual instrument becomes very simple. 
It is only necessary to provide for mounting the four fixed elements, the series elements of the 
T, in a base to which are attached the two sets of multiple terminals. 

As actually constructed, the instrument takes the form shown in figures 3, 5, and 6. The 
four series elements are mounted in holes bored parallel to the four sides of a 2 in. by 6 in. brass 
block. The holes are extended through the edges of the block to provide access to the com- 
ponents from one end, and to form the outer conductor of a short section of 3/4 in. coaxial line 
on the other. Holes bored through diagonally opposite corners from face to face provide the 
outer conductors for two more short sections of line. When center conductors are added, six 
short sections of coaxial line, whose center conductors meet in a point, are formed at two 
diagonally opposite corners of the block. Two of these line sections connect to the series 
components, and the other four are equipped with NBS-Woods [6] precision coaxial connectors 
to provide the bridge terminals. The junction point of the six line sections forms the “bridge 
corner. 

Capacitors C; and C, are coaxial capacitors with nominal values of 10 pF. C, is a parallel 
plate capacitor whose nominal value is 1 pF. RR, is a commercial metal film resistor whose 
value is selected to give the proper value of P at the frequency at which the instrument is to be 
used. 

Since it is necessary to adjust P, one or more of the series elements must be adjustable. 


C, is provided with two adjustments. Each plate has one half its surface cut about \%: in. 


lower than the other half, so that rotating one plate with respect to the other causes a change 
of about 10 percent in its capacitance. 


Final adjustment of P is accomplished by means of a 
trimming screw which extends into the space between the two plates. Gross changes in P can 
be accomplished by replacing the resistor R, with another of different value. 
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Figure 4. Equivalent circuit used in calculating 
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FIGURE } Cross section through the 
showing construction of the instrument. 


The basic instrument showing connector Figure 6. The basic instrument with standards 
attached for measuring a 200-ohm resistor. 


FIGURE 5. 
details and C, trimmer 


4. Experimental Results 


‘ l . : . , : 
P was adjusted to equal — at #=10' radians/sec, by measuring a 1000 Q resistor using a 


10 to 110 pF variable standard capacitor. This resistor, as were those used later, is a high 
quality commercial metal-film resistor mounted in a shielding can which is equipped with a 
NBS-Woods Precision Coaxial Connector. The three standard capacitors used are equipped 
with similar connectors. 
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When P had been adjusted to the point that the capacitance change on one terminal was 
within a few parts in 10° of that on the other, over an identical capacitance range, P was con- 
sidered to be adjusted to within a few parts in 10°. The adjustment was checked by measuring 
a 2000 © resistor using the same capacitor, a 20,000 © resistor using a 7 to 17 pF capacitor, 
and 100, 200, and 500 © resistors using a 100 to 1100 pF capacitor. The results are summarized 
in table 1. Since the error in resetting the capacitors is a few parts in 10° of their maximum 
capacitance, the adjustment of P is probably correct to within a few parts in 10°. 

With the values obtained for AC;, or AC2’, in checking the adjustment of P, and knowing 
that P is 10-7 + a few parts in 10°, it is possible to calculate the effective conductance of each 
of the resistors at 10° radians/sec. The results of the calculations appear in table 2. The 
measurements of C, obtained in checking the adjustment of P may be used to calculate a 
value for K. The results of the calculations appear in table 3. 


TABLE 1. Verification of the adjustment of P | TABLE 2. Comparison of measured resistance to d-c 


resistance 
% Diff 


AC? 
(ce, 


rermi 


pt 
999. 40 t pF Mho 
999. 37 1000, 44 1.00044 & 10 5 99. 971 
499. 71 4.9971 K 10 » , 200, 155 
199, 84 1. 9984 & 10 F 500, 428 
100, 031 1. 00031 XK 10 ¢ 999, 8453 
192. 92 50. 029 5.0029 & 10-4 ( 1999. 52 


192. 92 


499. 00 
498. 04 
5, 0180 5.0180 &K 10 19998, 4 
100. 050 
400.051 TABLE 3. Calculation of K jJrom me asured capacitance 
50. 206 
50. 208 
9994 


9991 
4. 9976 


The admittance measured by the instrument is the effective admittance at the “bridge 
corners.”’ In order to calculate Y, accurately, it is necessary to know the effect of the residual 
impedances between the terminal and the “bridge corner.”’ The equivalent circuit of figure 
4 is a reasonable representation of the residual impedances associated with the terminal. 

When the instrument is nulled initially with the terminal open, the effective admittance at 
the “bridge corner” is given by 


y G,1+R,G,) +e°R,Ci+ 3.o[C,—L,G—w’L,C7 
(1+ RG ,—oL 1)? T w(RC, T a0)" 


For the purpose of estimating the magnitudes of the residual impedances, the terminal can be 
considered to be 1% in. of 500, %in., brass air line. Since the series resistance and inductance 
in similar sections of line have been determined to be 0.8 milliohm/in. and 4.4><107>*° H/in. at 
10’ radians/sec, 2, is about 1.1 milliohms, and L, about 6.210-° H.C, is about 2.5 pF, and 
G, is less than 1 umho. Using these values for the residual impedances, the equation can be 
reduced to 

Y,.=G,4+joC,, at w=10' radians/sec. 


When an unknown admittance, Y,=G@,+jwC,, is connected to the terminal, the effective 


admittance at the “bridge corner’ is given by 


y GAG )N+ RAG +G) +? RAC+C,)?+jolC + C,—L (G,+-G,)?—o'L,(C+C,)?) 
1+ hA(G.4-G6,)—e7L (C,+-C,) P+e[ RCC, +L (6.4 4,)P 
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For G,<107? and C,<107"', this reduces to 


Y) =G,+6,4+jo(C;+C,—L.@) 
and 
aY,=Y—Y,=6,+ju(0,—L.@). 


The error in the measured conductance is insignificant, and the measured capacitance is 
in error by the term (—L,G,). Before the reactance associated with a nearly pure conductance 
‘an be accurately measured, it is necessary to evaluate this inductance, as well as the series 
resistance associated with the capacitor used to measure the conductance. It should be 
mentioned that these residual impedances cause no error in evaluating the constants of the 
instrument if care is taken to make all the terminals as nearly identical as possible. 


5. Discussion 


The measurements of P, which are summarized in table 1, indicate that it is possible to 
adjust P to within a few parts in 10° of 1/w. The limit on the accuracy of adjustment of P is 
the limit of stability and resettability of the standard capacitors used in making the adjustment. 

Comparison of resistance measured at 10’ radians/sec with that measured at zero fre- 
quency (table 2) indicates that the instrument is capable of measuring resistance to within a 
few hundredths of a percent. The very close agreement between rf and d-c values must be 
considered fortuitous at present because changes in capacitance are not known to much better 
than +0.05 percent. 

The values calculated for A (table 3) have a scatter of less than 1 percent of A. This 
indicates that it will be possible to make quite good measurements of the reactance associated 
with a nearly pure conductance. Apparently, the accuracy of such measurements will be 
limited by the accuracy to which the series resistance associated with the standard capacitors 
and the series inductance associated with the terminal can be determined. 

While the instrument has been proved only at a frequency of 10’ radians/sec, there is no 
apparent reason why the same principle can not be used at much higher, or at much lower, 
frequencies. The main deficiency of the instrument lies in the range of conductance which 
can be measured accurately. At present the range is limited, by the relationship G,=wAC, to 
two or three decades of conductance by the range of capacitance changes which is easily ob- 
tainable. Since the magnitude of the conductance measurable by a given capacitor depends 
on the frequency, the range of conductances measurable by the instrument is frequency de- 
pendent. The deficiency need not be too serious if it proves feasible to use the instrument to 
measure conductances for use as standards when calibrating other instruments capable of 
covering a greater conductance range. 

The instrument has demonstrated its ability to make conductance measurements limited 
in accuracy only by the accuracy to which changes in capacitance are known. This being the 
case, it is possible to immediately realize an order of magnitude increase in the accuracy 
of those conductance measurements which are within the range of the instrument. Of greater 
significance is the fact that the instrument demonstrates a principle which makes it possible 
to refer resistance measurements directly to standards of capacitance, thus placing the meast re- 
ment of resistance and reactance on a common basis. 


6. Appendix 


The null equations for the circuit of figure 2 are 


Y Y l > "i A) , sy ’ , ’ var y Y y ] 
GG; ( G, T Gstp )—w*(¢ 1G3+C3G,)(Cy+C;+ Cs) — w° CC ( G,- Gst+p ) 
6 t6 


G i'G ae ; l ’ ’ ve , ’ 
+3 (G,+Gs+-4;)- -a (p+ G,)(C— p+ OrtCy)— CC (G,4- G4 
6 6 


Wks 
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and 


] 


R, —w*( 10C3(C4 C;+Cs5) 


G,Gs(C.+-Cs+-C.) + (CG +CrQ) (G++ 


G, {~~ 1 ee C. HAS set aa ' 1 oa 
TR. ( ( srt 2-T Cs) +( Rt 1.) (G,+-G4,4+-G;)—w* 4 6 ( peered 27T C; )=0. (A2) 


A measurement can be made by connecting an unknown admittance, Y,=G@,+-jo(,, 
either to terminal-pair A, or to terminal-pair B. When Y, is connected to terminal-pair A, 
C, must be replaced by C3+-C,, G, by G:+G., G; by Gs, and C; by Cin the null equations. 
Similar changes must be made when the unknown admittance is connected to terminal-pair B. 
If AC,=Ci—C,, AG,=G:—G,, AC;=C;—C;, and AG;=G:—G;, and subscripts A and B 
denote the terminal-pair to which the unknown admittance is connected, subtracting the in- 
itial null equations from the final ones yields the four expressions 


2G sCe MO. Team Bi nr 
Coa | 14" = |--[1+45 |ac.+|* 7 C10; | a. 
C, C, wl, C, 


G y Y Y R, y Y 
+. id ah R,¢ | (AG, T G, )— ai ( ( iGs T ( "2G, JAC ‘c 
W 


on GF 4 


Cy w* Rol r¢ 3_ RG; a wo Rel 1C3__ ReGiGs AC. + 14 BeGsCe AG. 
C, ‘b ‘OF C's ( 


4 


Ry 
C, 


Tera wo RCC 5G, Gs G 
a. iz “ } BC s_ LieGh 73 AC; 4 ROa—— (AC.4 C,) 
Cs CU, bP C, 


GG 
- (¢ G34 CG, ) (AG: t G.)+( a ww Rf 'b ) AC AS 


4 


R '  , RG, 
pes -£ (C1G; 1 CsG,)4G.—( +s 
Cs ‘OF 


and 


WC, BGG, Gy, _ berngaing Ws. 
Gp [= ( ; -- ( | -—( 1 ') ac i (¢ 1G; r¢ 3G) (AC st Cz) 
4 


wf 4 4 ( 4 


3G: om GF OF ‘ G; : . 
+| A > _ - |: t . —— Ref | AG). (A6) 
wl, C, wC'; 


If ©, C3, and C, are good three-terminal capacitors, their losses will be due entirely to 
series resistance, in which case, 
G,=07Rh, Ci, G;~0R3C5, 
and 
G, = wh, ri 

If C; and C; are good two-terminal capacitors, changes in their conductance will be due entirely 
to series resistance, and 

AG. = wR2(C2—C?)=0?R,(Co+C3)AC2, 

AG; =~ wR; (C,+C.) AC, 


In the instrument described, 


C,~C;~10 pF 
C,=1 pF 
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C,<1000 pF 
C,<1000 pF 
R,~10002 at 10° radians/sec 
and R 


?, ~ 1002 at 10° radians/sec. 


It is not unreasonable to assume Ff, =~ Rh, = R, = Rh,~R;< 0.12 and 


a 


g~s 


rhe null equations can then be reduced by using the following approximations 


URC, YRCER Cs 
+ fig ‘|-[1+" R, ah ( |; 1t+wR,R.C,C;~1, 
C, C 


ee ~~ C. | ~— C: Sf? RRC, 1 "a _ Foil - 
wl" . ( 


4 


w Rel nCs _& . G; = be C3 RG Gs . 2 Rf id 3, 
( "4 w*( s Sie ( r 
G tC 
| Ss = — RP, 1 ay =-—R,¢ = Ry 1 — RC 
|G w Ref 1 |= 6 | Rel | . w| R46 on Rf 
4 


[¢ G3 T ( 3G, | ~~ [wR ¢ at 4 wR, a 3] ~~ wl, 1C: (B,C, 
Using expressions (a), (b), (d), and (f), eq (A3) is reduced to 


BOs e ROC 
— ; a = AG + (Ry 1—R,6 5) (AGo- é.,.)—— ee (R, 


4 


C,:+ RB3C;)AC;. 
4 
Since 


AG; = w R;(¢ ‘st Cz)AC 59 


AG, = wR2(¢ 'a-t ( < )AC 2 


ep “ils aan ’ 
wR, _ 7 AC 5 = Gea 


4 


Cra ~~ AC [1 w Ro(C, T C2) (RC, —R;¢ wI—-[RiC, T ?3¢ 3— EC, T RC, T BAC, T CVG x. A7) 


which is eq (16). 


Using (a), (ce), (e), and (f), eq (A4) is reduced to 
| 


, ; : . AC.( R,¢ 
1 (Ry ¢ i+ R;C;)(AG;+G_,,) + = ( 
( R, 
Since 
AC, 


0.0,> 
2, 
6 C, 


= Grr, 


wr | R,¢ 7 R3¢ 3) R;(¢ 
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this can be further reduced to 


( 7B == — AC st [RC 1 + R¢ 3 Ry ‘at Rol 6 Rs, ( 2 T C3) 1G, 8, 
which is eq (17). 
Using (a), (c), (e), and (f), eq (A5) is reduced to 


—w’ R,( i¢ 3 
C" 


SAC, —w*(R,C,— RyC,)(AC2+C,a) 


{ R,¢ 1 T R3( 3)AG;— AG». 


Since 
AG; = w R;(C;+ C3) AC; 
and 
w*( RC, + R303) R3(C;+ C3): 
this reduces to 
Y r Ref ‘ ( 3 y 9 y y Y \ 
Gui*— 7 — AC;—w*(R,C,— RC.) (402+ Cra) —AG?2 
which is eq (14). 


Using expressions (a), (b), (d), and (f), eq (A6) is reduced to 


Grp= 0 —w*(R,C,+ R;C;)(AC;+ Cx) —AG;+ 
R, ca 3 
( 4 


(R,C,— ReCe) AG 
CLC 
R, 7 3 


4 


AGi = wo Ro( C2 + C5) AC 2, 


9 T Cs) Ry "1; — PC 5) |—w*( BC, T R;C;) (AC; T Crz)—AG;, (A10) 


eq (15). Equations (A7), (A8), (A9), and (A10) should be accurate to within a few parts 
in 10° at frequencies up to 10° radians/sec for the instrument described. 
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Exact Inductance Equations for Rectangular Con- 
ductors With Applications to More Compli- 
cated Geometries 


Cletus Hoer* and Carl Love* 
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Exact equations are given for the calculation of the self-inductance of rectangular con- 
ductors and of the mutual inductance between combinations of parallel filaments, thin tapes 
and rectangular conductors. A general procedure is also given for calculating the self- 
inductance of complicated geometries by dividing the geometry into simple elements whose 
inductances can be calculated. This general procedure is valid for conductors having 
nonuniform, as well as uniform current densities. 


1. Introduction 


It is the purpose of this paper to give a number of exact equations for the self- and mutual 
inductance of rectangular conductors, and also to show how these equations may be combined 
to obtain exact inductance solutions for more complicated geometries. These equations were 
derived during a recent investigation of thin tapes as possible high frequency inductance 
standards. 

A number of approximate equations exist for calculating the inductance of rectangular 
conductors [1—5].'. The usual method is to calculate the mutual inductance between two fila- 
ments spaced a distance apart equal to the geometric-mean-distance [6] of the conductor or 
conductors. This method assumes that the length is much greater than the other dimensions. 
However, when the length is no larger than ten times the next largest dimension, the error from 
this method may be as large as several percent. It is for those cases where a high degree of 
accuracy is needed that the equations in this paper have been derived. 


2. Mutual Inductance Calculations 


In general, the mutual inductance between two conductors will be a function of the current 
distribution in each conductor. A conductor having a constant cross-sectional area? along 
its length may be thought of as a bundle of parallel filaments, each having a cross-sectional 
area (A and carrying a current Jd. The current density, J, is assumed to be constant along 
the length of each filament but to vary from filament to filament. The mutual inductance, 
M, between two conductors having constant cross-sectional areas A; and «A, and carrying 
currents J; and J; may be derived from energy considerations and is 


| . . 
M AF | 1; a MJ Jd Ad Ag, (1) 


where ./,,. is the mutual inductance between a filament carrying a current J,d.\, in the first con- 
ductor and a filament carrying a current J.d:A, in the second conductor. Since it is assumed 
that the current is constant along the length of each filament, the mutual inductance between 


*Radio Standards Engineering Division, NBS Boulder Laboratories, Boulder, Colo, 
Figures in brackets indicate the literature references at the end of this paper. 
L.e., tne cross-section dimensions do not change. 
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any two of the filaments is given by Neumann’s formula, 


(2) 


M,.=0.001 dl,-dly 
Ji Sle r 


The distance r is between two elements of length dl, and d/, which are along two filaments 
having total length /, and /,, respectively. For dimensions in centimeters, 1. is in micro- 
henrys. Since inductance cannot be calculated or specified except for a closed circuit, the 
integrals of (2) are usually over a closed path. However, the self- and mutual inductance of 
parts of a closed circuit (such as those to follow) may be calculated, provided that the total 
inductance of the circuit includes the contribution from each part. 

If the current density is constant throughout each conductor, (1) reduces to 


1 . | . » 
M- A,\A 5 | Madd, (3) 


which is independent of the current in either conductor and a function of only their dimensions. 
Using (2) and (3), the mutual inductance between any two conductors having constant cross- 
sectional areas and uniform (but not necessarily equal) current densities may be evaluated. 

Results obtained from (2) and (3) are quite useful even though the current densities are 
not uniform. In many cases, the conductors may be subdivided into elements small enough 
so that each element may be considered as having a uniform current density. The mutual 
inductances between the smaller elements are then properly summed to get the total mutual 
inductance. Results from (2) and (3) are also useful when the cross-sectional area of a con- 
ductor is not uniform along its length. In this case the conductor is subdivided into sections, 
each of which has a constant cross-sectional area along its length. The mutual inductances 
between sections are then properly summed to get the total mutual inductance. The technique 
of subdividing conductors into smaller elements whose inductance can be calculated is developed 
at the end of this paper. 


2.1. Mutual Inductance Between Parallel Filaments 


The first step in calculating the mutual inductance between parallel tapes or bars is to 
calculate the mutual inductance between parallel filaments using Neumann’s formula. The 
mutual inductance, /,, between two parallel filaments of length /,; and /, spaced in any relative 
position is 


M,=0.001 | z In (2+y2?+ 2) — ¥2+?? 


where 


Figure 1. Two parallel filaments whose mutual 
inductance is given by eq (4). 








The dimensions are defined in figure 1. If the left end of l, is to the left of the y-axis, the value 
of 1; will be negative. The value of M, is in microhenrys for dimensions in centimeters. This 
one equation covers all possible positions of two parallel filaments. It is given in its expanded 
form for special cases elsewhere in the literature [1, 7]. 


2.2. Mutual Inductance Between a Thin Tape and a Filament 


The mutual inductance, M,,, between a filament and a thin tape parallel to the filament 
such as shown in figure 2 is obtained by putting M, into (3) and integrating over the width of 
the thin tape. 


M,—' [“at,ar, 
i=q JM 


where p?=P?+ (H—7z)? in M,. 

Since z in M, is not a function of x (that is, the limits of integration s,, 82, s3, and s, are 
independent of x), it is unnecessary to expand M,, and the two terms inside of the brackets 
of (4) can be integrated as they are. This integration gives 


0.001] | 22—P? a > 
M,, ; A In (a+y2?+P?+27) + gz In (e+ y2?+P?4 
i" ( 


where 


(x) | (2) =>5355 (—1)***f (qi, &). 
" t , 1 : 
The positive directions for £, P, and J; are taken as the positive directions of the z, y, and 
> axes, respectively. If the back end of the filament is in any octant other than that shown 
in figure 2, some or all of the values of F, P, and J; will be negative. In this equation and in 


those following, the principal value of the inverse tangent is to be used. 


2.3. Mutual Inductance Between Two Thin Tapes 


To obtain the mutual inductance, /,, between two parallel thin tapes of zero thickness, 
such as shown in figure 3, M, is integrated over all filaments in both tapes. 


+y 





+2Z 


Figure 2. A filament parallel to a thin tape whose FiaurRE 3. Two parallel thin tapes whose mutual 
mutual inductance is given by eq (6). inductance is given by eq (8). 





1 °7 E+d fa 
M=—{|{ | Madr, 
; ad , E J0 — ‘i 
where p? +-(r4,—a2,)? in M,. This integration gives 


0.001 zr 2? P? 
V, 2 (2+ yz" t ——r In (4+42°+ P?4 
ad ; 2 


Lz 
Pyx’ +P? 


c? — Pz Tan 


where 


ss 4 4 
)=>0>0(—1)'***f (qi, s)- 


1k=1 


The distance F is measured from the yz, plane to the left edge of the second tape. If 
the back left corner of the second tape is in any octant other than that shown in figure 3, some 
or all of the values of &, P, and J; will be negative. 

The mutual inductance, 7,,, between two thin tapes whose axes are parallel but whose 
widths are perpendicular to one another such as shown in figure 4 is 


l 
a, 


M,, 


> Pic fa 1 ® P+ ¢ 
M,drdy M,,dy, 
= de Cc Jp 


where P in M,, has been replaced by the variable y. This integration gives 


0.001 a ; ; 
M, ( at )y In (2+y2?+y'4 


ac “ 


where 


i ~ +j+k+1¢ 
24 24 (—1) I (Gis Vis 8 
84 - 


rz plane to the bottom of the second tape. 


Figure 4. Two thin tapes whose aris are parallel 
but whose widths are pe rpendicular. 


The mutual is given by eq (10). 





2.4. Mutual Inductance Between a Bar and a Filament 


The mutual inductance, ./,,, between a rectangular bar and a filament parallel to the bar 
such as shown in figure 5 is obtained by integrating ./,, over all thin tapes in the bar. Equation 
(3) gives 
7p 


| 
b Jp 


Pe eee | | | " Mardy 
ab Jo Jo 


MAY, (11) 


where 7 in /,, is now replaced by the variable Y=P—y. The integration of (11) is identical 
to the integration of (9) except for the limits of integration and the constant. 


0.001 
ab 


M,, 


where f(z, y, 2) is the expression within the inner brackets of (10). The distance P in (12) is 
measured from the rz plane to the filament. 


2.5. Mutual Inductance Between Rectangular Bars 


To obtain the mutual inductance, .7,, between two rectangular bars, V/, is integrated over 
all thin tapes in both bars. 


ae |" | Mudrndys, 
be re 0 


where ? in M, is replaced by the variable (y.—y,). Since the x and 2 limits of integration on 
M, are independent of y, the four terms inside the brackets of (S) may be integrated as they 
are. After a lengthy integration, (13) yields an exact expression for the mutual inductance 
between two parallel rectangular bars spaced in any relative position. 


~4 


0.001 yet yk 2h , 
a abed ( 4 aaa a In ( 


1 ( ry =. 


Figure 5. <A filament parallel to a rectangular bar 
whose mutual inductance is given by eq (12). 








Figure 6. Two parallel rectangular bars whose 
mutual inductance is given by eq (14). 


a 
where St(z, y, 2) (x) Ty) (2) =>5 DS DS (Dt Ff, rj, 8). 


The dimensions are defined in figure 6. Either bar may be placed so that its back lower left 
corner is at the origin. The second bar may fall in any of the octants. If the back lower left 
corner of the second bar falls in any octant other than that shown in figure 6, some or all of the 
values of /, P, and /; will be negative. The distance # is measured from the yz plane to the 
left side of the second bar; P is measured from the xz plane to the bottom of the second bar, 
and /; is measured from the ry plane to the back of the second bar. 

The mutual inductance equations (4) through (14) have intentionally been left in their 
rather complicated form because it is in this form that one can most efficiently program them 
on a computer. For most cases, considerable effort is required to evaluate these expressions 
without the aid of a computer. If the dimensions vary greatly in magnitude, the positive 
terms may very nearly cancel with the negative terms, and care must be taken to calculate 
each term to an accuracy sufficient to assure the required accuracy in the result. 

It should also be noted that if either x or y or z approaches zero, all inverse tangents in (6) 
through (14) go to zero.’ If any two of the variables x, y, and 2 approach zero, all terms in 
(14) go to zero except the square root term. 


3. Self-Inductance Calculations 


The self-inductance of a conductor is a special case of the mutual inductance between two 
conductors. We can think of the self-inductance of a conductor as the mutual inductance 
between two identical conductors which coincide with each other. The discussion on mutual 
inductance calculations may then be applied to self-inductance calculations. The self-induct- 
ance, L, of a conductor having a constant cross-sectional area, A, and carrying a current J is 
given by (1) where the two integrations are now over the same conductor. 


=}, f f. Med Jd Ad Ay 


The mutual inductance, 1/,., between two filaments in the conductor is given by (2) if the current 
JdA in each filament is constant along the length of the filament. 
If the current density is constant throughout the conductor, (16) reduces to 


b=, [ f MudAd Ay, 


which is a function only of the dimensions of the conductor. 


3 In eqs (6 





For the self-inductance, L,, of a thin tape of width “a” and length “J, (17) gives 


L, " ‘ | : M dz,dro, 


where p?=(2.—2,)? in .M,;. This integration gives 


0.002 [+ ~vl?+-a? : ee a+ y¥a?+[ 
L,.=—S> | 0% In . — (P+a*)?? + 3a ln . + /°+a? |- (18) 
Od” a 


This result may also be obtained from (8), the expression for the mutual between two thin 
tapes, by setting d=a and E=P=1;=0. 

The self-inductance, L,, of a rectangular bar of sides ‘‘a’”’ and “6” and length ‘J’ is most 
readily obtained from (14), the expression for the mutual inductance between rectangular bars, 
by setting d=a, c=b, and E=>P=/1,=0. These substitutions give 


0.001 
a*b? 


where /(x, y, Z) is the expression within the inner brackets of (14). It can be shown that f(z, y, 2) 
is an even function of z, of y, and of z This result allows one to use the absolute value of 
the limits of integration, thereby reducing (19) to 


0.008 


a’h? 


L, 


where 


mn Js ; 
f(x, y, 2) 2) Kw) (z)=>5 >> D> (—1)'t*** fai, 7), 8) 
1 1 k=1 


12 r2 $2 


For bars with a thickness-to-width ratio less than 0.1, the self-inductance L, 


may be 
calculated from the empirical equation 


LO", 21) 


where L, is calculated from (18). The values of A given in table 1 are sufficiently accurate 
to calculate Z, to an uncertainty of 1 ppm for the range covered by the table. These values 
were obtained by forcing (21) to give results equal to that obtained from (20). 


TABLE 1. Values of K, uh/cem 


0. 002 
. 005 
O10 
O15 
020 
030 
040 
050 
060 
O70 
USO 
O90 
100 


OSY 2. OSY 
OS2 2. 082 
0712 OF17 
OO17 O624 
0529 2 O53S8 
0364 0376 
0210 0226 
0066 | 2. 0OS5 
9929 Q9OAD 
9797 S21 
9071 697 
9548 


¢ 9577 
9430 M461 
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FIGURE 7. Two arbitrary circuit elements car? ying 
currents i, and ip. 


4. Inductance of Complicated Geometries 


The inductance of complicated geometries or of conductors having nonuniform current 
densities often can be calculated from a consideration of the magnetic energy of the conductor. 
Consider the two circuit elements shown in figure 7. The instantaneous energy, w, stored in 
the magnetic field created by the sinusodial currents 7, and 7, may be written 


w= bL ii + gLets+ Mitr, (22) 


where L, and L, are the self-inductances of elements “1’’ and ‘2’, respectively, and 1), is the 
mutual inductance between these two elements. In (22) the currents are the instantaneous 
values 


inn SiN (wt +¢,), 


where i, is the maximum value of 7, and ¢; is the phase angle. 
[t is often more convenient to use the complex notation 


T,=a,+)b,, [2=a2+jbo. 


Using this notation, the average magnetic energy W’ stored in the system of figure 7 may be 
written 


W bLl, dT; if s Lele: Ts T Mol, - Is. (23) 


It is assumed that the currents J, and /; have the same frequency but not necessarily the same 
phase angle. The pseudo vector notation /,- 7, has been used to take care of the difference in 
phase angle. 


Linton 


I, - 2=4,d2 4 bib, — COS (91 — 2). 


Equation (23) for the average magnetic energy of two circuit elements may also be written in 
summation form as 
~ La 
(2: 
W=;5> 25 Mil / (24) 
where ./;,=self-inductance of element 7. 
For a circuit having n elements, the average magnetic energy of the total circuit is 
1S 


W.== > DM yl, - I. 


The equivalent inductance, L,, of a closed cireuit having n elements is that inductance 
which will give the same average energy, W’,, when the same total current, J7, is flowing through 
the circuit. 


W, tLwlr: Ie. (26) 
Solving (25) and (26) for the equivalent self-inductance gives 


| 5. 
leelr i T= 


M,,1, + I;. 


i 
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FIGURE 8. A conductor of complicated cross-sectional FIGURE 9. 
area divided into elements whose inductances can be 
calculated. 


-An equivalent circuit of the conductor 
shown in figure 8. 


This equation can be used to calculate not only the equivalent self-inductance of a closed 
circuit but also the equivalent self-inductance of part of a circuit if there is negligible inter- 
action between the part under consideration and the rest of the circuit. The self-inductance 
of conductors having complicated geometries, therefore, may be calculated from (27). The 
procedure is to subdivide the conductors into simple geometries whose inductances can be 
calculated, and then sum these values with (27). 

For example, let us calculate the equivalent self-inductance of the conductor shown in 
figure 8, assuming that there is negligible interaction between it and the rest of the circuit. 
There are no formulas for the self-inductance of such a geometry, so divide the conductor into 
elements whose inductances can be calculated. Since the self- and mutual inductance of 
parallel rectangular bars in any position can be calculated, divide the conductor into n=7 
rectangular bars. The conductor now looks like the circuit in figure 9 with n=7 elements in 
parallel. The voltage drop across each element can be written 


n 
e=Tent jo 3 Mud, (28) 
i=1 


where r, and M,, are the resistance and self-inductance of element *, and /,; is the mutual 
between elements & and 7. Since the desired self-inductance is independent of the value of e, 
set e=1+ 0 for simplicity. Equation (28) gives n equations which may be solved for the 
n unknown currents J; through J,. With the currents known, the equivalent inductance can 
be calculated from (27); 


where, for this example, 


At low frequencies (28) becomes 
in which case (29) reduces to 
(30) 
where F? is the total resistance of the conductor. Since the conductor is assumed to have a 


constant cross-sectional area along its length, each element will have a common length, /. 
The total resistance, 2, may be calculated from 





where A; is the area of element & which has a resistivity p,. Equation (30) shows that at low 
frequency the inductance of a conductor is not only a function of its dimensions but also of 
the resistivities of its different elements. When the resistivities of all elements are the same, 
however, (30) reduces to 


Sy ee 


Pig F l t= 
which is a function of only the dimensions. If the 3 sections of the conductor in figure 8 all 
have the same resistivity, the low frequency inductance may be calculated directly from (32) 
where the total area A is 


The procedure outlined in this example is valid for a conductor of any arbitrary constant 
cross-sectional area, provided of course that the conductor can be subdivided into n elements 
whose inductances can be calculated, and that there is negligible interaction between this 
conductor and the rest of the circuit. 

If the interaction is not negligible between the part being calculated and the rest of the 
circuit, the whole circuit must be considered. If there is interaction between the conductor 
of the previous example and, say, m other conductors in the circuit, the sum from i=1 to n in 
(28) becomes i=1 to n+m. There would then be n+™m currents to be evaluated, after which 
the equivalent inductance of the conductor could be calculated from (29), i and 7 going from 
1 to n as before. 

* If exact values of inductance are not required, the equivalent self-inductance often may be 
calculated easier using geometric-mean-distances [1-7] (abbreviated g.m.d.). It should be 
remembered, however, that even exact equations for the g.m.d. gives approximate values of 
inductance because all g.m.d. formulas assume that lengths are infinite. 


5. Conclusions 


Exact equations have been given for the mutual inductance between the following parallel 
conductors: a filament and a thin tape, a filament and a rectangular bar, two thin tapes, and 
two rectangular bars. There is no restriction on the size or spacing of the conductors. Exact 
equations have also been given for the self-inductance of a thin tape and of a rectangular bar. 
Since these equations are quite long, they have been put in a form most suited for programming 
on high speed computers. 

These equations for simple rectangular conductors may be used to calculate the self- 
inductance of any other more complicated geometry if that geometry can be divided into simple 
rectangular elements. A general procedure for calculating the self-inductance of complicated 
geometries has been given. The mutual inductance between complicated geometries may also 
be calculated using a procedure similar to that used for self-inductance. 

In deriving the equations given in this paper, the assumption was made that the current 
density was uniform throughout the conductor (or throughout each element in the case of com- 


plicated geometries). This restriction is not a severe one however, since in many cases the 
conductor may be divided into elements small enough so that the current density in each 
element may be considered uniform. Applications of these equations to the calculation of 
high frequency inductances where the current density is not uniform are given in a forthcoming 
paper by Brooke, Hoer, and Love [8]. 
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The volume common to two cylinders of radii r;<r. with axes intersecting at angle 
8 is found to be rjv(k)/sin 8, where k=r,/r, and v(k) may be evaluated (1) as the hypergeo- 


metric series 
x fs 1 
4, 45 23k*)=4e D(C 2) am 
292, a, sy 
. n/\ 1 


n 


2) as the combination of complete elliptic integrals (8/3)[(1+ k?) E(k) 1—k?)K(k)] or 
3) as the cumulative integral 


k 
8 | kE(k)dk. 


A table of v(k) to 8 decimals over the range 0<k(0.01) <1.00, including 62, modified second 
central differences, is presented. This volume integral was useful in interpreting a gamma- 
ray albedo experiment involving a collimated source and a collimated detector, and may also 
be applicable to crossed-beam experiments. Two series useful for k close to unity are pro- 
vided, one of which involves differencing against the series 


16/3=4" D7 (7 )(* 1): 


1. Introduction | for this region, derived from the right-angle elliptic- 


| integral solution [4, 5, 6], is more complicated but 
| also more rapidly convergent. 


In crossed-beam experiments [1]! using the high- 
intensity accelerators now becoming available, the 
“ceometrical target,’’? or volume common to the two 
colliding beams, is a useful parameter for interpreting 2. Volume Integral 
the measured data. An evaluation of this volume in The 
terms of an infinite series was recently exhumed for : 
possible application to an x-ray free-air ionization 
chamber having a cylindrical sensitive volume inter- 
sected by a pencil of x rays [2]. This evaluation 
had been used in the analysis of a gamma-ray beam 
back-scattering experiment [3] for making a theoreti- 
cal estimate of the single-scattered component of 
the radiation ‘‘seen’”’ by a collimated detector. 

Evaluations of the volume common to two circular 
cylinders of unequal radii with axes intersecting 
at right angles [4, 5, 6], and of equal radii with axes 
intersecting at an arbitrary angle |7], have frequently sinB J, 
been offered as calculus textbook exercises. How- 
ever, a combined treatment does not seem to appear 
in the technical literature in a form cin for 3. Common Volume When r,=r; 
easy application to practical problems. The follow- 
ing results provide formulas, a table, and a graph For equal cylinder radii r;=r2=r, 
for such applications. (1b) reduces to the familiar result [7] 

The series used in [3] is here corrected, expressed 
in terms of binomial coefficients, and identified as a 7 "a, (ss 
hypergeometric series. For nearly equal cylinder V(r, B) sin B (r—2")dx 
radii, convergence can be accelerated by use of the ti 
difference-series technique [8]. An alternative series 


integral for the common volume of two 
cylinders of radii 7;Sr2 with axes intersecting at 
angle 8 (see fig. 1) is found as follows. The cross 
section parallel to the cylinder axes, at a distance 
z from them, is a parallelogram of height 2(73—22)'/? 
and base 2(77—<2?)'//sin 8B. Hence the volume inte- 
eral is 


(1b) 


the integral in 


167° 


1 Figures in brackets indicate the literature references at the end of this paper. | 3 sing 
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FIGURE 1.* Three-view sketch of the common volume of cylinde rs 
with*radii r, and 1, axially intersecting at angle B. 


area of the shaded parallelogram parallel 
eft view cc 


to the plane of the axes in the lower 
ymprises the integrand in eq (la) and is integrated over the range 
r, shown in the other two views 


4. Common Volume When r,<r, 
4.1. Series Solution 
The factor (73—2°)' in the integral in (1b) may 
be expanded as a power series |9, p. 2, eq 5.3] 


in z/ro, since xSr,Sro. The volume integral then 
becomes 


2, B) al ( 


sin 6B *< 0 \tt 


ae edn a 
;) Gi-2) "dr 


which may be integrated term by term. 
The resulting series solution 1s 


sin B . l 


This same result may be obtained by casting the 


| integral in (1b) in the form 


V(r, 72, B)=(r3/sin B)4 ( = ) 
“ 


| 0 t)'? (1-3 dt, (6) 
0 Io 


where t=22/r7, which is now recognizable as the 
integral representation of the hypergeometric series 
[10] 


Vir, r2, 8) =(ri/sin B)2n (7) 2F, (—3, 5 2: 
re e 2 
identifiable with (5). 

Since (5) is somewhat slowly convergent when 
r,;&ro, under some circumstances it may be advan- 
tageous to difference this series against a 
series (16) discussed in the appendix, giving 


Arr 4 
sin B 


The convergence rate of the series-term in (8) is 
not improved over that of (5). However, for 
r,&r, this sum is small compared to the constant 
term 4/37, identifiable with the equal-radii solution 
(3), hence resulting in higher precision of V(r, 
ro, B) for the same number of terms. An alterna- 
tive series solution for this region is given at the 
end of the following section. 


1 /x- 


Vir, ro, B) 


4.2. Elliptic Integral Solution 


An alternative solution of the integral in (1b) may 
be obtained as a combination of complete elliptic 
integrals [11] of the first and second kinds, A(k) 
and E(k). Applying formulas (219.11) and (361.03) 
from Byrd and Friedman [12] the result is found to 
be? 


; Sr ee ase 

1, 72, B)==—— Lk?) E(k) —-(—Fk 9 
Viri, r2, B) 3sing kh?) E(k (1—k*)K(k)|] (9) 
where 


\/ le. 


Except for the angle factor 1/sin 8 this result is 
the standard textbook solution [5, 6] for cylinders 
intersecting at right angles. Also, the formulation 
in (9) is related to the indefinite integral [12, eq 
(611.01)] 


. 
( kE(k)dk : [(14+h)E(k)—U—Fk)RK(k)). 


2 This integral is part of the ‘‘G factor’’ used for interpreting gas scattering ex- 
periments in which a circular-aperture detector views a gas target transversed 
by a cylindrical beam. In this context this elliptic integral solution has been 
given by E. A. Silverstein, Nucl. Instr. and Meth. 4, 53 (1959) and by D. F 
Herring and K. W. Jones, Nucl. Instr. and Meth. 30, 88 (1964 
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A partial check on (9) is obtained by expanding 
E(k) and K(k) as power series in k according to 
reference [9], eqs (773.1) and (774.1). Combining | Tasve 1. Values of v(k), defined in eq (12a-e), over the range 
like powers of & and substituting back r,/r2 for k the 0 Sk(0.01) $1.00, valid to the 8D given 
resulting series is identical with (5). An additional | Modified second central differences 32 are provided for interpolation u 
check is provided by the identity of eqs (9) and (3) nuiiary tables 
in the limit as k->1. Also, eq (9) can be obtained 
from (7) by use of the tenth Gauss recursion formula 
on page 9 of reference [10]. 

For k close to unity, a series which converges more 
rapidly than (8) may now be derived by substituting ices | oe 
in (9) the series in equations (773.3) and (774.3) in nica co 
reference [9] for A(k) and £(k). The first few terms 05 | 0.0157 0305 
of this series are 


in which values of v(k) are interpolated from table 1. 


0.00 0. OOOO 0000 +2 


06 0226 0928 
0307 6874 
0401 


y > . 16 d ? O508 
V(r), 72, 8) =(73/sin B) -— , 
») 


0627 


O759 


4 1\,4 . VP 2 | 0.008 


~ 1400 
Z. __ (In : 


2°-4°-6°-8 


3 & 1603 
—_ ) k’ 7 1809 

, * yor 
kv AC ment 


2204 


2500 


32.5. 4 967\,,10 
—4: rp. (In5—= ys sca (11) 1 0.2755 


3022 
3301 
: 3592 
where k’”°=1—k?=1—(7/r.)?._ Using only the terms 

given in (11) the sum for 7;/r2=0.90, without the 

factor (73/sin8), gives 4.49991482. as compared 

with the exact value of 4.49991288 .. ., and the 
convergence improves as 7)/fe goes toward unity. 

The series in (11) may be obtained in general form, 

if desired, by use of the appropriate transformation 

[13] on the hypergeometric series given in (7). 


5. Numerical Results 


In table 1, the dimensionless factor 


rik) 


sin B 


I's 


Vir, Te, B), (12a) 


(12b) 


(12c) 


s{ kE(k)dk. (12d) 


J 0 


Modified second central differences 62, are included 
for interpolation by Everett’s formula 
s(1 +h) B(k)—(—k) K(k), (12e) , , 
= | v(k,; +pAk) =(1—p)v(ki) + potkiss) + E2(p)&,, ; 
where k=r,/rz is tabulated to 8 decimal places for ; 
0<k(0.01)<1.00, computed using (12e) and K(k) + F2(p)8m, vit +++, (14) 
and E(k) from [11]. Hence, for many practical 
applications, the common volume of two cylinders | where pis the interpolation fraction of the interval of 
with radii 7,;<r. and axes intersecting at angle 6 | tabulation Ak, and 6, ;, &, ,., are the modified second 
may be computed as differences at the tabular points 7 and 7+1 and were 
evaluated from the second and fourth differences 
7 ey according to 
V(r, %2, 8) == v(k) (13) | 
/ sin B 





5? =6’—0.1846'. 
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Everett’s coefficients £,(p) and F2(p) are available | | An additional 1/z-series, which also does not appear 
in standard tables [14, 15] and are identical with | in standard compilations of series [17], can be ob- 
Lagrangian interpolation coefficients A*, and A} | tained by combining (16) with a series discussed by 
[16, table 25.1]. Bromwich [18], [17, eq 274] 
The general behavior of v(k) is shown by the curve | 

in figure 2. Values of x(k) can be taken directly ¢=o/ity 
from this curve for use in rough calculations where - > 
only two- or three-figure accuracy is required. 


to form 


0 


= — 
heal 
n=1 


n—1 


Using the addition theorem 


( a )+( nes l 


the result is 


Figure 2. Graph of v(k) showing the general behavior of the 


, The author thanks E. Hayward for suggesting 
function and suitable for rough calculations. = . we wat Mines Ben — 5 
aii the problem, and P. Lamperti and L. W. B. Jolley 
for their stimulating interest in the cylinder inter- 
section and 1/m series, respectively. The author is 
. . also indebted to F. W. J. Olver, I. A. Stegun, L. F 

. Appendix. Two Series for 1/x — i ee ee ee es 
6. App / Epstein, and A. Fletcher for suggestions and com- 
ments, and especially to Mrs. Ruth Capuano for 


The series evaluation (5) for V(7r,,72,8) contains a : 
Sy : ~ | computing table 1. 


factor of z and the formula (3) for V(r,8) does not. 
Thus, for r;,=r2=r, the right-hand side of (5) can 
be equated to the right-hand side of (3) to form aie 
16r 4nr } A 
= ; “2 be (15) [1] A. Sehoch, A discussion of colliding beam techniques, 
3sing sing Nuclear Instr. and Methods 11, 40 (1961). 
[2] P. Lamperti, private communication. 
from which [3] E. Hayward and J. H. Hubbell, The backscattering of the 
» Co® gamma rays from infinite media, J. of Appl. 
Phys. 25, 506 (1954). 
|4] W. E. Byerly, Elements of the Integral Caleulus (Ginn 
and Co., Boston, 1902), p. 281, ex. (5a), which is 
incorrect. 
[5] H. Hancock, Elliptic Integrals p. 89, ex. (6) (John Wiley 
(16) & Sons, New York, N.Y., 1917). 
[6] P. Franklin, Methods of Advanced Calculus, p. 299, ex. 77 
(McGraw-Hill, New York and London, 1944). 
aca : = “o~ [7] See, for example, J. Edwards, A Treatise on the Integral 
Phis series can now be used to form the difference- Calculus. I., p. 793, ex. (23) (MacMillan, London, 
series 1n eq (5). 1921). 
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This paper presents a mathematical analysis of two-dimensional steady-state heat 
conduction in a solid of rectangular cross section, two parallel surfaces of which are exposed 
to separate temperature environments with respective surface heat transfer coefficients 
For the other two parallel surfaces, the temperature environment and surface heat transfer 
coefficients are assumed to vary as a function of position on these surfaces. A particular 
application of this analysis has been to determine the temperature distribution in a trans- 
verse cross section of exposed concrete columns. 


1. Introduction 


Modern architectural treatment of exterior columns of high-rise buildings involves the 
use of columns exposed to temperature extremes on the outside face and relatively constant 
temperatures on the interior.’ The temperature extremes give rise to dimensional changes in 
the column which can best be determined from temperature distribution in the column. In 
connection with a request from the Federal Housing Administration, an investigation was 
made for determining the temperature distributions in an exposed exterior column. 

This paper presents a mathematical analysis of two-dimensional steady-state heat con- 
duction in an exposed exterior column which may be represented by a solid of rectangular 
cross section, two parallel surfaces of which are exposed to separate temperature environments 
with respective surface heat transfer coefficients. For the other two parallel surfaces, the 
temperature environment and surface heat transfer coefficient are assumed to vary as a function 
of transverse position on these surfaces. 

An analysis for general application is presented, along with the statement and numerical 
results for four specific problems. 


2. Analysis for General Problem 


Because the heat-transfer problem in a rectangular solid may have a more general applica- 
tion than that cited in the previous section, the following analysis is based on boundary con- 
ditions, arbitrary in nature, which may be adapted to specific applications. 
the assumed boundary conditions for heat transfer in an exposed exterior column. 


Section 3 presents 


For steady-state, two-dimensional heat conduction in a rectangular homogeneous solid 
with no flow of heat in the direction of the third dimension, it is assumed that the rate of heat 
flow at each of the four boundary surfaces is proportional to the temperature difference be- 
tween the temperature of the surface and that of the ambient adjacent to the surface. The 
proportionality factor may be prescribed as a function of the distance along the surface. A 
convenient placement of coordinate axes for the type of problem to be considered is as shown 
in figure 1, where the boundary conditions at the four surfaces, r= +1, and y= +a, are shown 
adjacent to the surfaces. 
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FicurE 1. Transverse section of rectangular solid, 
showing boundary conditions on x 1 and 
Vv ta, 


The temperature potentials V and V(r) are the temperature differences between the 
temperatures of the ambient adjacent to various surface positions and that of the ambient 
adjacent to r=l; the temperature potential v is the temperature difference between a tempera- 
ture in the solid and that of the ambient adjacent to z=/. The proportionality constants 
H, and H, represent the coefficients of heat transfer between the surface and the ambient at 
the surfaces s=—l, and r=l, respectively. For the purposes of this paper the variable heat 
transfer coefficient /7(r) and temperature potential V(x) are considered to be piece-wise con- 
tinuous on —l<xr<l, where H(r)=H,(r) and V(r)=V;,(r2) are continuous on z;_;<r<z; and 
i=1,2,3, ...m;2=—l, and z,=. 

The partial differential equation for steady-state heat conduction in a homogeneous solid 
in a rectangular coordinate system is 

Ov Orv 
Oat Oy - (1) 


A solution of (1) satisfying the boundary conditions at r=/ and r= — (fig. 1) is 


vr hi Ayhol 1—2z/l) >_> A,R(B,2 l) cosh B,Y l 
Vi hy thet+2hyhe ’ = cosh £6,a/l 


° . °,° i dv ; 
where due to the similar conditions on y= +4, it Is assumed that - 0 on the center line y=0. 
( 
Symbols used in (2) are: 
K=thermal conductivity of the solid 
AIK 
ALK 


coefficients to be determined from boundary conditions on y=a 
B, cos B,(1+2/l)-4 h, sin B,.1+2 1) 


where 8, are positive roots of (6?—h,h2) sin 28B—B(hi +h.) cos 28=0. 
Substitution of boundary conditions on y=a gives m equations 


a lH (zx) lLH(x) FT V(r) hi t+hho(i—z/l) 
Ss (2B x ‘ S u 
A, R(8,. | 8. tanh B,a/l+ m4 K h, thet 2hihs ] 


Multiplying each equation by F(8,2/l) and integrating from z=—l to r=1, yields 


Sz 


>>} BA, tanh B,a/l R(B,2/l) R(B,x drt 73 As | 
n=1 J —! i=l e 


H(xr) R(B,27/l) R(B,x Daz | 


Zi-1 


i Hx) R(Gx nx" 3 ) hy t+hoh,(1—ar 2) de. (4) 
1 


/ m °r 
ee | h, T ho T Qhyhe 


{- 
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Using the orthogonal and other related properties of ?(8,2/l) as derived in the appendix, (4) 
may be put in the form of simultaneous equations 


B,Ant+ 2) Ch, rAb Di, 
k=1 
kAn 
where 


7 
1 


l a Zi 
B,=8, tanh B,a/l ( (8.2 drt p> [ H (x) R?(B,,2/l)dxr 
sat Jt) 


1 


i 


H(z) R(B,2/l) R(B,.2/l)dr 


| st Vir) Athlke—sez 2] 
> ce : rag ay. 
D, K 7 1 I [ V hy, ho ZQhihe H(x)I (Bua l) ly 


Of particular interest is the mean temperature in the solid, v,,, given by 


l . [* 
Sal f . i} y dude 


h,A+he) l > A, tanh £6,a//[B, sin 28,+h,(1—cos 28.) |. 
hi thet+2hhe  2aG B*. 


3. Temperature Distribution in Exposed Concrete Columns 


An interesting specific case for solution of (2) is to consider a transverse section of a rec- 
tangular concrete column. Herea portion of the boundary surface is exposed to outside weather 
conditions or design temperature, another portion of the surface is exposed to the ambient 
temperature maintained within the building, and another portion of the surface is in contact 
with the walls of the building. A cross section of the column, and an abutting wall in the region 
b<xr<d, is shown in the lower portion of figure 2. The r=/ face of the column is assumed to be 


VARIATION OF H (x) AND V (x) 
FOR -KPeH (x) [ v-V (x) ] ON y=a 


Ficgure 2. Transverse section of rectangular concrete 
column with abutting wall separating the inside and 
outside environment. 


Outside face is located at r=1. 





the outside face exposed to the outside design temperature which is assumed to be the zero or 
datum plane temperature from which all other temperature potentials are reckoned; the face 
r / is assumed to be exposed to the indoor temperature potential V. 

The upper portion of figure 2 shows the assumed variation of //(r) and V(x) along the 
surfaces y=+a. H(x) is assumed constant at 77, in the region —/<xr <b; decreases linearly 
from H, to zero in going from 6 to ¢; increases linearly from zero to //, in going from ¢ to d; 
and is constant at HZ, in the region d<r</. V(r) is assumed constant at V from —/ to ¢; 
and is zero from ¢ to/. In the region of the abutting wall, d<cr<d, the heat transfer between 
the wall and the column was assumed to be dependent on the inside and outside temperatures 
in their region of influence and the surface coefficients decreasing linearly to zero at r=c. 


de 
—K dy 


H,(< 


Expressions for D,, B,, and C,,, necessary for substitution in (5) are given in the appendix 
(A13, Al4, A15). The solution of (5) for A, and its substitution in (2) gives the temperature 
distribution for values of x and y. 

Numerical solutions were obtained, using a digital computer, for a 36 in. 14 in. concrete 
column for four positions of the abutting wall, based on a 100 deg F temperature difference 
between the inside and outside. Numerical values used, and assumed as constants, were 


Thermal conductivity, A=1.0 Btu/hr ft deg F 
Inside heat transfer coefficient, 77,—0.5 Btu/hr ft? deg F 
Outside heat transfer coefficient, 7,=6.0 Btu/hr ft? deg F 
Width of column, 2a=14 in.=7/6 ft 
Thickness of column, 2/=36 in.=3 ft 
Thickness of wall, d—b=8 in. =2/3 ft 
Also it was assumed that d—c=c—6. 


The inside surfaces of the column were assumed to be separated from the inside of the 
building by suitable interior finish, giving a nominal value for //, of 0.5 Btu/hr ft? deg F. The 
value H,=6.0 Btu/hr ft? deg F was based on an outside wind velocity of 15 mph. 


MEAN TEMP 
33.4° 


INSIDE AMBIENT 
TEMP 


Figure 3. Temperature potential distributions in 36-in. by 14-in. concrete columns for four positions of the 


abutting walls, based on a 109-deg temperature difference between the inside and outside and the constants 
cited in the text. 
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Figure 3 shows temperature potential isotherms plotted within the column cross section 
for four positions of the abutting wall. Mean temperatures were computed from (6). Fig- 
ure 4 shows mean temperatures over the width of the columns plotted against distance meas- 
ured from the inside surface of the column. Positions of the abutting wall are shown for 
clarity. 

Temperatures in degrees Fahrenheit can be expressed in terms of the temperature potentials, 
v, of figure 3 by means of the relation 


t, °9F=t,+(v/100) (t, 


where ¢, is the inside air temperature, and f¢, is the outside air temperature, in degrees 
Fahrenheit. 


4. Discussion 


The analysis of section 2 assumes similar boundary conditions on the surfaces, y= +a. 


For the consideration of a more general case where these conditions are not similar, it is con- 
venient to make a translation of the z-axis of figure 1 to the lower side and rewrite the sum- 


mation term in (2) 
») 
E cosh B,y/1+E, cosh B, (= y | R(8,2/L) 


cosh B,2a/l 


Substitution of boundary conditions at y=0 and y=2a will give two sets of simultaneous 
equations similar to (5), both containing the coefficients A, and &,. Solution for A, and 
FE, may be obtained by iteration; for example, let #,,=0 in one set and solve for A,, substi- 
tute A, in the second set and solve for /,, which is substituted in the first set. This process 
may be continued until A, and £, are invariant. 


Figure 4. Mean temperatures over the width of the 
columns against distance measured from the inside 
surface of the columns. 


Positions of the abutting walls are shown 
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The assumption was made that a sufficiently good approximation to the solution of the 
infinite set of equations (5) in the infinitely many unknowns A, could be obtained by solving 
the first forty equations for A,, n=1,2,...,40. Solution of the first twenty equations for 
A,, n=1, 2..., 20 gave results which differed from that for forty A, by less than 0.1 °F. It 
therefore appeared that convergence of the approximating sequence of A,’s to the true solution 
of (5) was sufficiently rapid to justify our assumption. 

A computer program in Fortran IV language is available, for the problem as set forth in 
section 3. Data input to this program allows any possible variation in the numerical values 
for the thermal conductivity, surface coefficients, and dimensions. 


4. Appendix. Orthogonal and Other Related Properties of /?(8,,7//) 
To satisfy the boundary condition in (2) at r=/, 8, must be the roots of 
(B?—hyhz) sin 28—B(h, +h.) cos 28=0. (Al) 


Without proof in this paper, the roots 8, will be assumed real and positive. Proof of the real 
nature of these roots is similar to that for roots of an expression given in Carslaw and Jaeger.’ 
It then becomes expedient to develop the orthogonal properties of 


R(6,xr/l)=8, cos B,(1+27/l) +h, sin B,11+27/l) (A2) 
where —/<.r<l. Several expressions useful for later developments are cited below: 


B?(hy +he)? 
(6?+-h7)(6°+h2) 


sin® 28 
d (B?—hyh2)? 
(6° +hi) (8? +h) 


cos- 28 


R?(Bx/1) a Fe cos 28(1+2/l) +48 sin 28(1+-2/l). 


R(B,2/l) R(B,2/L) Ea(2)- La(2) 


) 


L(x) =(8,8;—h?) cos ¢@(1+2/l)+hid sin (14+ 2/1) 
L2(r)=(8,8,+h?) cosy (1+2/l)—h,wW siny (14+ 2/1) 


o= 6,-4 Bx, y= B,»—B, 
M,(z) [ R(Gu2 l)dz 5 (8, sin B,(1+2/l)—h, cos B,(14 


N, (2) | R?(B,2/l)dx mens hi) eo sin 26,(1+2 ye cos 28,(1+2/l) 


Py (x) L,(x)dr 5 ((Bu8:—Ai) sin (1+ 2/1) —hy@ cos ¢(1+-2/L)] 


Qo(x) L.(xr)dr i [(8,8,+h7) sin ¥(1+2/l)+hy cos ¥(1+2/1)|. 


Carslaw, H. S., and Jaeger, J. C., Conduction of Heat in Solids, 2d ed., p. 114 (Oxford Unive rsity Press, London, 1959). 
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By use of (A3) and (A4), it ean be shown that Pol) + QoQ = Po(—L) + Qo(—L) =0; Le., the 


integral R(8,2/1) R(B,.2/l)dr=0, showing that R(8,7/l) and 2(8,7/l) are orthogonal to each 


other in the interval —/<r</. Also, using (A3) and (A4), 


N,(1) : (82-+h?)(B2+h?+he) 


- (87 +h) 


l 
—No(—D=5 (Brthith). 


Assuming /7,(2) and V(r) to be polynomials in sz, then the recurrence 


M,(z) zr? R(Br/l)dr=2M,_\(r)— | M,_\(2)de All) 


and similar relationships for N,(2), P?,(2) and Q,(2) may be established for determining B,,; 
D,, and (,,, in (5). Assuming a resulting polynomial of the form S+ 7r+ Uz? in the solution 
for )),, the indefinite integral becomes 

l 


g2) Mo(x)+ 7 (T+2U'r) R(Byx/l). (A12) 


. 9 2 ’ 
S+7Trt+U 2’) R(B,2/l)dz ( S+T7r- ue—2tt 


The solution for D, for the specific case illustrated in figure 2 takes the form 


LTiaett x?) R(B,2/l)dx 


p= f 


hytho+2hihe 4. 


where 


1=0 
S,)/(e—b) -T,/(e—b 
T,/(d—e) 


=0. 


17 
1 U 
S,—cT,)/(d—c) U 
7 


The function in the region —/<(r</ is a linear, piecewise continuous function for which 
the integrals go to zero at the integration limits 7 l and r=l, and elsewhere at the other 
limits only the coefficients of /?/82 in (A12) need be considered, 


D, liyhs (2/‘U.M,(b)+(U;—U.) M,(e)—U,M,(d) 


Br (hy+he+2hyh2) 
+ (7,+ T,—2U 2b) R(B,,6/l) —(T,— 7T3—2U 3d) R(B,,d/l) 
T;,—T.+: 
By a similar analysis 
(8, tanh B,a/1+h2)(B2 +h?) (82 +h3+h2)+(8, tanh B,a/l+h,)(8,+h3) (8, +h? +h,) 
2(82-+-h?) 


hl 


: ] hol »2 »2 1 »2 > »2 | 
1433 E = R?(8,d/l) —R?(B,e/l) ee R*(g,,¢/1) —R*(8,,b/L) 


ad 


&; |. ft?) hile : 6) he(d +) 
4 l l 


k hol L,(d)—L,(ce) L;(d)—L.(c) hil L,(c)—L,(b) , L2(e)—L2(b) e 
2(d at o T y Has | ros 7 y } (A15) 
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More wavelengths from thorium lamps, W. F. 
R. W. Stanley, J. Res. NBS 69A (Phys. 
(Mar.—Apr. 1965), 70 cents. 

For the purpose of supplementing or superseding the spectro- 
scopic secondary standards of wavelength derived since 1910 
from an iron arc, Meggers and Stanley, in 1958, reported the 
first interferometric determinations of wavelengths emitted 
by a thorium lamp. Those determinations were restricted to 
222 intense radiations of thorium with wavelengths ranging 
from 3288.7356 A to 6991.5839 A. Now the same, and addi- 
tional, interference spectrograms have been measured to 
provide improved wavelengths for 511 radiations ranging in 
wavelength from 3269.6089 A to 7020.504 A. The present 
list includes many lines of lower intensity than those pre- 
viously published and fills most of the large intervals in our 
first report. The accuracy in relative wavelength values of 
163 classified thorium lines is tested by the combination 
principle which indicates that the average error is less than 
1 part in 20 to 40 million. Similar measurements of wave- 
lengths emitted by iron-halide lamps have errors that are 
3 to 5 times greater. 


Splitting of a set of equivalent sites in centrosymmetric 
space groups into subsets under homogeneous stress, J. B. 
Wachtman, Jr., and H. 8S. Peiser, J. Res. NBS 69A (Phys. 
and Chem.), No. 2 (Mar.—Apr. 1965), 70 cents 

It is assumed that the symmetry elements possessed by a 
homogeneously stressed crystal will be those common to the 
crystal and to the macroscopic state of stress. Application 
of stress either leaves a space group unaltered or lowers it to a 
subgroup. Such lowering can always be considered to take 
place in successive steps each of which leaves no group 
intermediate between the starting group and the subgroup. 
Each such step can be accomplished by a uniaxial stress; 
for the centrosymmetric space groups, all but one of the 
symmetry reductions consisting of two or more successive 
steps can also be accomplished by a uniaxial stress. <A set 
of points all of which are equivalent in the unstrained erystal 
frequently splits into two or more subsets under stress. For 
each space group all possibilities are taken into account by 
considering the behavior of the general position because the 
behavior of each special position can easily be derived by 
specializing the description of the general position. The 
splitting of the general position under homogeneous stress is 
tabulated for all possibilities for each of the crystallographic 
space groups possessing centers of symmetry. 


Self and mutual admittances of waveguides radiating into 
plasma layers, J. Galejs, Radio Science J. Res. NBS/USNC 
URST 69D, No. 2, 179-189 (Feb. 1965) $1.00. 

The self and mutual admittance of waveguide backed slots 
radiating into a plasma (or dielectric) layer has been for- 
mulated in a laterally unbounded geometry. The ad- 
mittance expression involves integrals that may be approxi- 
mated by summations for numerical calculations. For a 
constant mesh size approximation to the integrals the sum- 
mations are the same as obtained in earlier work, where the 
plasma layers are considered to be within a wide waveguide. 
The accuracy of the solution is improved by decreasing 
the mesh size used for the numerical integration, which is 
equivalent to increasing the size of the large guide in the 
waveguide model. 

Mutual admittance calculations between two-parallel slots 
show that the presence of a plasma layer decreases the mutual 
admittance relative to its free space value. 





Electroacoustic waves excited by a space vehicle in ionized 
atmosphere and its effect on radar return, K.-M. Chen, 
Radio Science J. Res. NBS/USNC-URSI 69D, No. 2, 
241 (Feb. 1965), $1.00. 
When a conducting 


~~ 
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space vehicle in the ionosphere is 
illuminated by an incident electromagnetic wave from the 
ground, a current and charge are induced on the vehicle. 
The induced current and charge in turn generate two scattered 
fields in the ionosphere; an electromagnetic (EM) and an 
electroacoustic (EA) wave. A part of the EA wave is con- 
verted to an EM wave across the discontinuity of the iono- 
sphere and reaches the ground. As the result, the Radar 
return of a space vehicle in the ionosphere is greatly enhanced. 
In this paper, the case of a conducting cylinder in a plasma 
and illuminated by an EM wave is considered. The induced 
current and charge on the cylinder are determined and the 
scattered EM and EA waves are calculated. The cor- 
responding radar cross sections due to EM and EA waves 
are defined. It is shown that the enhancement of the radar 
return due to an EA wave bears some resemblances to the 
large outbursts of the reflected HF signals from the satellites 
observed by Kraus. 


Insulated and loaded loop antenna immersed in a conducting 
medium, R. H. Williams, Radio Science J. Res. NBS/USN( 
URSI 69D, No. 2, 287-289 (Feb. 1965), $1.00. 

Previous analyses of the driving-point impedance of a small 
loop in an insulating radome immersed in a conducting me- 
dium are extended to include the effect of a permeable core 
loading the loop. It is shown that the effective area of the 
loop is increased by a gain factor previously derived for the 
loaded loop in free space. The effects of the conducting 
medium in this factor are negligible. 


Electromagnetic waves along an infinitely long and thin 
conducting wire in a magneto-ionic medium, Y. Mushiake, 
Radio Science J. Res. NBS/USNC-U RSI 69D, No. 503- 
510 (Apr. 1965), $1.00. 

Klectromagnetic fields of an infinitely long conducting cylinder 
in a magneto-ionic medium with axial static magnetic field 
are theoretically analyzed. Expressions of the electromag- 
netic fields and the dispersion formulas are obtained. Ex- 
plicit approximate expressions of the relative propagation 
constants for the case of an extremely thin wire are derived 
and their dependence on the plasma and the cyclotron 
frequencies is discussed. The behavior of the electromagnetic 
fields around the wire is also discussed. 


ts 


Excitation of acoustic waves in plasmas, 
Radio Science J. Re s. NBS USNC 
616, (Apr. 1965), $1.00. 

This paper reports early results in an experimental study to 
determine the effects of perturbing weakly-ionized gaseous 
plasmas with acoustic waves emanating from transducers 
that operate in the audio and ultrasonic frequency ranges. 
Compact loudspeakers incorporated into cylindrical plasma 
discharge tubes provide the source of acoustic signals. Mod- 
ulation of electromagnetic waves by the acoustically disturbed 
plasma is measured in a unique rectangular cavity whose out- 
put is detected and fed to a wave analyzer. Resultant 
wave-analyzer responses indicate that the collision frequency 
is modulated in addition to the plasma frequeney, and 
suggest that the variations in both are proportional to the 
magnitude of loudspeaker diaphragm deflection, as predicted 
by simple acoustic theory. Knowing that two forms of 
modulation occur simultaneously, a calibration scheme was 
developed to separate the variation Af, in the plasma fre- 
quency from the total response. Measurements to date 
indicate that AN,./N.=WAN/N (where N. and N are electron 


W. A. Saxton, 
URSI 69D, No. 4, 609 
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and neutral-molecule densities, respectively) for slightly- 
ionized gases which are subjected to low-frequency pressure 
variations in the order of 10-5 mm Hg, and that such varia- 
tions produce plasma-frequeney perturbations of 0— 140 ke/s 
at plasma frequencies up to 800 Me/s. 


Test of the constancy of the velocity of electromagnetic 
radiation in high vacuum, P. Beckmann and P. Mandies, 
Radio Science J. Res. NBS/USNC-U RSI 69D, No. 623 
628 (Apr. 1965 $1.00. 

It is pointed out that Einstein’s postulate of the constant 
velocity of light is verified only indirectly by elementary- 
particle experiments leaning more or less heavily on present 
electromagnetic theory, the latter being verified only for low 
velocities. Direct experiments can be explained by the bal- 
listie theory of light if transparent media, such as gases, re- 
radiate as a secondary source. A direct experiment with 
coherent light reflected from a moving mirror was performed 
in vacuum better than 1076 torr 
the constant velocity of light. 


Its result is consistent with 


Signal degeneration in laser beams propagated through a 
turbulent atmosphere, P. Beckmann, Radio Science J. Res. 
VBS/USNC-URSI 69D, No. 4, 629-640 (Apr. 1965), $1.00. 
The statistical distributions of the angle of arrival, the spot 
location, the cross section, the amplitude, the carrier phase, 
and the modulation phase of a laser beam traversing an 
anisotropically turbulent atmosphere are derived in terms of 
the space correlation function of the atmospheric index of 
refraction and the windspeed ‘he limitations imposed by 
the turbulent atmosphere and the loss of coherence on the 
depth and bandwidth of the modulation, on the length of the 
path and on the aperture of the receiving apparatus are 
analyzed. Experiments to obtain numerical parameters and 
to check functional dependencies are proposed. 


Consolidated Supplement to Building Materials and Struc- 
tures Report 144 (Supersedes Supplements 1 and 2 of BMS 
Report 144). Sound insulation of wall, floor, and door con- 
structions, R. D. Berendt and G. E. Winzer, NBS Mono. 7? 
Vov. 30, 1964 40 cents. 
The data obtained at the National 
the sound insulating properties of 
summarized. The results of the 
to BMS Report 144 


results obtained through January 


Bureau of Standards on 
building structures are 
two previous Supplements 
included, together with later 
1964. Single figure ratings, 


1955) are 


STC and INR, for airborne sound transmission and impact 


sound transmission, respectively, as well as the octave fre- 
quency band spectra of impact noise, are included as addi- 
tional information. A brief description of the sound-measur- 
ing techniques is given. 


Investigation of the hydraulics of horizontal drains in plumb- 
ing systems, R. S. W vly, VBS Mono. 86 (Dec. 18, 1964 
OO cents, 
tesults are reported from an investigation of the hydraulics 
of flow in experimental apparatus simulating nominally 
horizontal simple and branching drains of plumbing systems. 
The data are correlated with limited findings in an earlier, 
unpublished NBS study the results of which have been 
utilized in current plumbing codes. The need for further re- 
search is pointed out, particularly in relation to hydraulic 
performance of drain systems as affected by steep slopes, 
drain storage volume, energy losses at stack bases, attenuation 
of water depths and discharge rates in long drains, and large 
drain diameters. 

Analysis vielded equations useful in estimating hydraulic 
capacities for surge flow or for surge flow superimposed on 
steady flow over a range of conditions. Capacity estimates 
for the large drains are substantially greater than formerly 
assumed. The experimental findings suggest that through 
careful design it may be possible to reduce the number of 
conventional vents extending above plumbing fixtures in 
small systems such as those in one- and two-family houses. 
The results, of interest to writers of plumbing codes and 
handbooks, have important applications in computation of 
optimum plumbing loads on building drains and building 
sewers, in terms of actual plumbing fixtures. 


Ellipsometry in the measurement of surfaces and thin films 
Symposium Proceedings Washington 1963, ed. KE. Passaglia, 
R. R. Stromberg, and J. Kruger, NBS Misc. Publ. 256 
(Sept. 15, 1964), $2.25. 

A Symposium on the Ellipsometer and its Use in the Measure- 
ment of Surfaces was held at the National Bureau of Standards 
in September 1963. This volume contains 19 of the papers 
included in the program, together with any discussions which 
followed oral presentation. Topics covered include historical 
review, theory, computational techniques, measurement 
techniques, and the use of ellipsometry in measuring metal 
surface oxide films and organic films. 


Standard Reference Materials: Metallographic characteri- 
zation of an NBS spectrometric low-alloy steel standard, 
R. E. Michaelis, H. Yakowitz, and G. A. Moore, NBS Misc. 
Publ. 260—3 (Oct. 22, 1964), 20 cents. 

The spectrometric standard steel designated NBS Low-Alloy 
Steel 461 was investigated by means of electron probe micro- 
analysis and quantitative metallographic techniques employing 
a digital computer. Electron probe microanalysis showed 
the steel to be homogeneous in nickel and iron at two to four 
microns of spatial resolution. The average of all determi- 
nations agreed with the certified values for these elements. 
Inclusions in the steel were identified, classified as to size 
and shape, and counted. Mean free path data on the inclu- 
sions were calculated. The ASTM ferrite grain size number 
was deduced as 13.5 for the steel in the unetehed condition. 
From the mean free paths in ferrite and pearlite, it was found 
that the steel is structurally homogeneous at a five micron 
level. It is concluded that the homogeneity level corresponds 
closely to the grain size of the material. It is further con- 
cluded that NBS—461 steel is sufficiently homogeneous that 
any present microanalytical technique can be carried out 
with little chance of inaccuracy due to inhomogeneity. 


Legibility of alphanumeric characters and other symbols, I. 
A permuted title index and bibliography, D. Y. Cornog, F. C. 
Rose, and J. L. Walkowiez, NBS Misc. Publ. 262-1 (Dec. 15, 
1964), 60 cents. 

This permuted title index and bibliography to the literature 
on the legibility of alphanumeric characters and other sym- 
bols includes 325 references Lu 
lap of terminology in the legibility literature, this report 
uses the term, “legibility’’ to include “‘legibility,” ‘‘read- 
ability,”’ “perceptibility,” ‘‘visibility’’ and any other closely 
related concepts. “Other symbols’? include such items as 
the arrows and other coded symbols used to present informa- 
tion on radar displays. The psychological literature on per- 
ception has been included only when it was closely involved 
with the specific problems of alphanumeric characters and 
other meaningful symbols. Studies concerned with en- 
vironmental variables, i.e. illumination and symbol-back- 
ground contrast, and the legibility of dials and scales have 
received little attention. An author index is ineluded. 


) 


to the confusion and over- 


Aluminum nails, CS263-64 (Nov. 1, 1964), 10 cents. 

This Commercial Standard covers the types, sizes, material, 
and principal dimensions of aluminum wire nails commonly 
used in the building industry. The finish, the approximate 
number of nails per pound, and tolerances for the principal 
dimensions are included. Provision is also made for the 
labeling or identification of nails complying with the standard. 
The standard does not cover specific applications or special 
purpose nails. 


The electrical properties of aluminum for cryogenic electro- 
magnets, R. J. Corruccini, NBS Tech. Note 218 (Aug. 30, 
1964), 30 cents. 

The published data for the ideal resistivity and the mag- 
netoresistivity of aluminum have been correlated. It is 
shown that both properties can be calculated for the limited 
ranges of temperature and purity that are of importance for 
aluminum cryogenic electromagnets from the residual resis- 
tivity ratio alone. Empirical functions are given for these 
properties, and sample calculations are given of the figure of 
merit for a particular aluminum solenoid relative to the same 
solenoid using water-cooled copper. 
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Radio path length stability of ground-to-ground microwave 
links, M. C. Thompson, Jr., and H. B. Janes, NBS Tech. 
Note 219 (Nov. 15, 1964), 45 cents. 

The Lower Atmosphere Physics Section of the National 
Bureau of Standards, Troposphere and Space Telecommunica- 
tions Division, conducted a series of eight field experiments 
in 1958 and 1959 to study the variation in the apparent path 
length of ground-to-ground microwave links (i.e., the varia- 
tions in phase of 9400 Me/s signals) and the corresponding 
variations in radio refractive index and other atmospheric 
parameters measured at the path terminals. This report 
summarizes the results obtained on paths in Colorado and 
Florida. It also includes an analysis of long-term variations 
in the apparent path length and associated variables recorded 
during each experiment, and a power spectrum analysis of 
short-term variations in apparent path length and refractivity. 


Some modifications in methods of calibration of universal 
ratio sets, ). Ramaley, NBS Tech. Note 220 (Aug. 30, 1964), 
25 cents. 

Universal Ratio Sets can be calibrated by a number of differ- 
ent methods. The well established methods are very briefly 
outlined and emphasis is placed on some more recent devel- 
opments. The choice of methods will depend upon available 
laboratory equipment and other considerations. 


Bibliography on the measurement of bulk resistivity of semi- 
conductor materials for electron devices, J. C. French, NBS 
Tech. Note 232 (Oct. 21, 1964), 60 cents. 

In support of a study of accurate measurements of the bulk 
resistivity of semiconductor materials such as germanium 
and silicon which are used in electron devices, a literature 
search has led to the development of a rather large collection 
of references. The file cards making up the bibliography are 
reproduced in this Technical Note for the assistance of others 
who are concerned with the improvement and standardization 
of these measurements. 


Survey of magnetic thin film materials, G. W. 
NBS Tech. Note 247 (Se pt. 25, 1964), 40 cents. 
This survey lists the materials reportedly made as a magnetic 
thin film, along with some of their properties and potential 
applications. Research activity using the less-frequently 
mentioned magnetic film materials is noted. About 200 
references are cited. 


teimherr, 


Research on crystal growth and characterization at the 
National Bureau of Standards January to June 1964, ed., 
H. F. MeMurdie, NBS Tech. Note 251 (Oct. 19, 1964), 46 
cents. 

The National Bureau of Standards with partial support from 
the Advanced Research Projects Agency of the Department 
of Defense is continuing a wide program of studies involving 
crystalline materials. These include investigation of methods 
and theory of growth, study of detection and effects of 
defects, determination of physical properties, refinement of 
chemical analysis, and determination of stability relations and 
atomic structure. The types of materials range from organic 
compounds, through metals, and inorganic salts to refractory 
oxides. This Technical Note, the fourth in the series, sum- 
marizes the progress of these various projects from January 
to June 1964, and lists the related publications and partici- 
pating scientists. 


Procedures for precise determination of thermal radiation 
properties November 1962 to October 1963, J. C. Richmond, 
D. P. DeWitt, and W. D. Hayes, Jr., NBS Tech. Note 252 
(Nov. 20, 1964), 40 cents. 

The preliminary design of an intergrating-sphere reflectom- 
eter, utilizing a helium-neon continuous-wave gas laser as 
the source, for measuring the reflectance of specimens at high 
temperature, was completed. Development work on an 
ellipsoidal mirror reflectometer for measuring spectral reflec- 
tance in the wavelength range of 2 to 15 microns of specimens 
at room temperature was continued. The study of equations 
relating spectral emissivity of metals to other properties was 
continued. Platinum-13°% rhodium and oxidized Inconel 
working standards of normal spectral emittance were cali- 





brated over the wavelength range of 1 to 15 microns at 
temperatures of 800, 1100 and 1300°K. Several modifications 
of the normal spectral emittance equipment were made to 
permit operation in the 15-35 micron range. 


Disclosures on various subjects; A plumbing vent manifold, 
a micro-adjuster, a glass joint or stopcock, and a miniature 
puller, ed., D. Robbins, NBS Tech. Note 253 (Nov. 30, 1964 
10 cents. 


This Note presents four devices embodying interesting and 
unusual solutions to problems prevalent in their respective 
arts. These devices, developed at the National Bureau of 
Standards and the U.S. Weather Bureau, comprise a vent 
manifold for a plumbing system; a micro-adjuster providing 
purely rectilinear motion; a miniature gear, flange and bearing 
puller, and a glass joint or stopcock preventing capillary 
action. 


On the formulation and numerical evaluation of a set of 
two-phase flow equations modeling the cooldown process, 
S. Jarvis, Jr. NBS Tech. Note 301 (Jan. 1965), 55 cents. 

A model of transient two-phase flow in a pipe is constructed 
in Eulerian coordinates assuming a single velocity but inde- 
pendent temperatures in the two phases. Experiments on 
the numerical integration of the system for cooldown prob- 
lem by both Lax and Courant-Isaacson-Rees methods indicate 
that a very fine spatial difference net must be used to com- 
pensate for the numerical diffusion essential to computational 
stability if a second surge is to be realized. 


National standard reference data system. 
tion, EK. Brady and M. B. Wallenstein, 
(Dec. 30, 1964), 16 cents. 

In June 1963 the National Bureau of Standards was assigned 
the responsibility by the President’s Office of Science and 
Technology of administering a government-wide program to 
promote and coordinate systematic data compilation and 
evaluation activities in all fields of the physical sciences. 
This government-wide program has been called the National 
Standard Reference Data System. The present report 
describes the background of this program and the mode of 
approach adopted by the management of the National 
Bureau of Standards in fulfilling its responsibility. 


Plan of opera- 
NSRDS-NBS1 


Weights and measures—a challenge to industry, M. W. 


Jensen, Proc. 50th Mid-Year Meeting Chemical Specialties 
Manufacturers Assoc., Chicago, Ill., May 15-20, 1964, pp. 
24-26 (Chem. Specialties Mfg. Assoc. Inc., New York, N.Y.). 
The responsibility for regulatory control over commercial 
weights and measures in the United States is left principally 
to the individual States and their political subdivisions. The 
National Bureau of Standards provides the States with 
calibration of the reference standards and various technical 
services pertaining to a weights and measures program. An 
important vehicle employed by the Bureau in this area is 
the National Conference on Weights and Measures. 
provides a national forum for discussion and a 
for developing model laws and regulations. 


which 
mechanism 


Computer interpretation of English text and picture patterns, 
R. A. Kirsch, EEE Trans. Electron. EC-13. 
No. 4, 363-376 (Aug. 1964). 

This paper considers a class of information sources consisting 
of text and pictures. The text is English language text 
appearing in scientific and technical documents. The picture 
sources are the largely schematic pictures that occur in the 
same class of documents. However, the discussion is 
broadened slightly to include other picture sources. For 
a tiny fragment of English, the paper shows how the syntactic 
structure of text may be described, and then goes on to suggest 
that a similar analysis may be performed on the class of 
pictures under study. The description of these two kinds 
of information with a single class of descriptive 
techniques is suggested as an alternative to the synthetie 
approach in which artificial languages are specified and then 
learned and used. The major reason discussed here for doing 
syntactical analysis of such sources is that several informa- 
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tion processing operations, amounting to the interpretation 
of the information sources looked upon as languages, can be 
done by the technique of syntax direction which uses the 
results of syntactic analysis to mediate subsequent processes 
for manipulating the information tokens. The paper con- 
cludes with an illustration of an algorithm for matching the 
sentences given by a simple grammar against the class of 
simple pictures which these sentences purport to describe. 
Methods of calibration of airglow photometers at Fritz 
Peak Observatory, L. L. Smith and R. B. Alexander, 1QSY 
Instruction Manual No. 6: Airglow, pt. II, pp. 25-36 (IQSY 
Secretariat, London, England, 1963). 

The problems of establishing low-intensity standard sources 
for calibrating airglow photometers are discussed; particular 
attention is paid to a radioactivated phosphor source which 
has desirable characteristics as a field standard. Attention 
is called to the necessity of adequate control of the trans- 
mission of any interference filters and the spectral response 
of any photomultipliers used in such photometers. 


Ozonation of asphalt flux, P. G. Campbell and J. R. Wright, 
I. & E. C. Prod. Res. and Devel. 3, 186 194 (Se pt. 1964 
The effects of ozone on the oxidation of asphalt flux were com- 
pared to those of air and of oxygen. Oxidative changes were 
measured by infrared spectroscopy, increases in softening 
point, and elemental analysis. 

In the oxidation of one flux of Southeastern U.S.A. origin, 
ozone enrichment of air reduced the time required to reach 
a given softening point to one-half that needed for air alone. 
Ozone-enriched oxygen reduced this time to one-sixth that 
needed for air. Reaction temperatures did not exceed 
158 °C, approximately one-half those used in commercial 
air-oxidation: gas flow rates were about one-fourth the com- 
mercial rate. Coating-grade asphalts produced under these 
mild conditions were more stable to photooxidation than was 
a commercial asphalt from the same flux. Similar results 
were obtained with a California-Coastal asphalt flux. 
Ozonation of a ferric chloride-catalyzed flux of California 
Coastal origin produced different results in that no oxidative 
changes could be detected in the residual hardened asphalt; 
the rate of hardening, however, was increased over that of 
other oxidants. Possible reactions explaining the role of 
FeCl; are presented. 


On the validity of the Einstein-Smoluchowski theory light 
scattering, R. Zwanzig, J. Am. Chem. Soc. 86, 
(Mar. 30, (1964). 

The molecular basis of the Einstein-Smoluchowski 
scopic theory of light scattering is investigated, using a 
generalization of Frohlich’s theory of dielectrics. The macro- 
scopic and molecular theories differ by an experimentally 
negligible quantity, associated with the electric field de- 
pendence of the dielectric constant. A slight depolarization, 
not accounted for by the macroscopic theory, is associated 
with the same quantity. Because of the method of deriva- 
tion, the results obtained here are valid only in the limit of 
long wavelength and low frequency. 
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Pulsed and CW sinusoidal voltage and current measure- 
ments, M. C. Selby, Proc. URSI XIVth General Assembly, 
Tokyo, Japan, Sept. 1963, XIII, 92-100 (1964). 

A Triannual progress review on the art of measuring single 
frequency and wideband voltages and currents in the USA is 
presented for the URSI. This is the first report on these 
quantities in the series of URSI, Commission I reports and 
covers the approximate period of 1960-1962. Present day 
significant applications, parameters and ranges of pulsed 
voltages and currents are briefly pointed out. Only the rela- 
tively more important publications are indicated and the 
more outstanding contributions commented on at some 
length. 


Thermodynamic data compilation and review at the National 
Bureau of Standards, G. T. Armstrong, Svensk Kem. Tidskr. 
76, No. 8, 419-427 (1964). 

Three examples of activities at the 
Standards involving compilation and 


National 
critical 


Bureau of 
review of 





thermodynamic data are described in some detail: (1) Se- 
lected values of chemical thermodynamic properties, (2) phase 
diagrams for ceramists, (3) thermodynamic properties of a 
selected group of light elements. The type of literature 
search, the scope of substances covered, the types of data 
gathered, the method of review and the manner of publication 
are described. A brief description is given of plans for the 
new Standard Reference Data Program at the National 
Bureau of Standards. 


Propagation of a reflected shock wave in a transverse mag- 
netic field, Y. Boublil and K. B. Earnshaw, Phys. Fluids 7; 
No. 6, 876-881 (June 1964). 
Shock waves are reflected from 
transverse magnetic field. The gas is highly ionized and 
therefore the lines of force can be considered ‘frozen in.” 
In such a case the magnetic field is compressed and the 
Rankine-Hugoniot relations, including magnetie pressure and 
magnetic energy as given by de Hoffmann and Teller, are used. 
The magnetic probes show an increase of the magnetic field 
behind the reflected shock, and the velocity of the reflected 
shock waves (measured from streak camera photographs) 
increases with increasing fields. For a constant incident 
velocity, the dependence of the reflected shock velocity on the 
magnetic field is given by a linear function: U,?=a+bB? 
where a and b are constant in a first approximation. This 
law has been verified by experimental results. 


a wall in the presence of a 


Refractivities of the rare earth oxides, E. M. Levin (Proc. 
Third Rare Earth Research Conf., Apr. 2-24, 1963), Book, 
Rare Earth Research, Chapter on Magnetic and Optical Studies, 
pp. 339-350 (Gordon and Breach, Publ., Inc., New 
York, N.Y., May 1964). 

Index of refraction and density data obtained for a number 
of rare earth borates (ABO,) and aluminates and gallates 
(A,B;O,2) permitted the calculation of the specific refractivi- 
ties (by the Gladstone and Dale formula) of the rare earth 
sesquioxides. This was done using the additive rule between 
the products of the specific refractivities and the weight frac- 
tions of the constituent oxides. The specific refraectivity 
values decreased with decreasing ionic radius of the cation, 
ranging from 0.148 for La,O, to 0.107 for Dy.O; to 0.090 for 
Lu,O;, and include the lowest values ever reported for oxides. 
Specific refractivities and molar refractivities (Lorentz- 
Lorentz formula) were calculated also from data on crystal 
fragments of the cubie rare earth oxides, obtained from two 
sources. Specific refractivities decreased from a value of 
0.120 for Dy,O; to 0.099 for Lu,O;, and were 10 to 13% 
higher than those calculated from the binary oxide com- 
pounds. Molar refractivities (Lorentz-Lorentz) of the pure 
oxides varied from 22.5 for Dy,Q, to 20.2 for Lu,O; and were 
in good agreement with reported data. 


Scie nce 


Analysis of a differential phase shifter, 1). A. 
IEEE Trans. Microwave Theory Tech. 
453-459 (July 1964). 

This paper presents the theory and analysis of a ganged pair 
of “line stretcher’? microwave phase shifters. The error 
analysis shows that some of the errors inherent in a single 
phase shifter of this type can be reduced through the use of a 
differential system; however, the magnitudes of other errors 
may more than offset the reduction. Graphical data are 
included to facilitate the rapid determination of the limit of 
error for any specified angle measurement. 


Ellerbruch, 
MTT-12, No. 4, 


Comparison calibration of inductive voltage dividers, R. V 
‘ A 
$< 


Lisle and T. L. Zapf, 1SA Trans. 3, No. 3, 238-2 
1964). 

Calibrated inductive voltage dividers serve as_ excellent 
standards of voltage ratio for use in the comparison calibra- 
tion of other a-e voltage dividers. Differences in phase 
angles in the dividers can be accommodated by the produc- 
tion of a voltage in the detector circuit in quadrature with 
the reference voltage. The paper discusses the comparison 
circuit, the handling of stray impedances, guarding, measure- 
ment procedures, quadrature voltage production, and errors. 
The comparison method described can be used in standardiz- 
ing laboratories for the calibration of inductive voltage 
dividers and other voltage ratio equipment. 
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Heater lead problem in calorimetry, D. C. Ginnings and E. D. 
West, Rev. Sci. Instr. 335, No. 8, 965-96? (Aug. 196 4). 

A study has been made of possible calorimetric errors resulting 
from the use of potential leads at the mid points of current lead 
segments between a calorimeter and its surrounding shield. 
With simple assumptions believed justified in most calorim- 
etry, there is derived a mathematical solution for this error 
in the case where the calorimeter and shield are at different 
temperatures and there are different thermal contacts between 
the current lead and the calorimeter and shield. It is con- 
cluded that when accuracies of better than 0.10 are desired, 
careful consideration should be given to heater lead design. 


Aspects of the National Bureau of Standards instrumentation 
literature reference file, J. F. Mayo-Wells (JSA 19th Annual 
Conf. and Exhibit, New York, N.Y., Oct. 12-15, 1964), ISA 
Preprint 20.6-5-64 (1964). 

The Instrumentation Literature Reference File 
associated Peek-a-Boo retrieval system of the National 
Bureau of Standards Instrumentation Division are briefly 
described. Some capabilities of multi-aspect indexing which 
extend beyond those of classical indexing are discussed, as are 
advantages of the so-called inverted file. For instrumentation 
information and related questions, the NBS Instrumentation 
File activity a means of access to the Bureau as a 
whole, which may in itself be regarded as a measurements 
information center of particular interest to other standards 
laboratories. Various NBS publications series are described. 
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Versatile Stark waveguide for microwave spectroscopy, 
D. R. Lide, Jr., Rev. Sct. Instr. 35, No. 9, 1226 (Sept. 1964 

A simple waveguide design for Stark-modulation microwave 
spectrometers is described. The guide provides a homo- 
geneous Stark field which is perpendicular to the microwave 
electric field. The design has been used successfully in 
spectrometers operating at high frequencies, high tempera- 
tures, and for other special purposes. 


Behavior of coaxial cable connectors for pulses with nano- 
second risetimes, R. B. Schwartz and A. C. B. Richardson, 
Nuclear Instr. Methods 29, 83-88 (North-Holland Publ. Co.., 
Amsterdam, The Netherlands, 1964) 

Reflections and risetime distortion 
mismatches were measured for several common eable- 
connector combinations, using a fast pulser and sampling 
oscilloscope (system risetime=0.4 nsec). It was found that 
excellent connectors are available for 1252 (RG 63/U) eable, 
as well as for 50Q cables. In addition, details are given for 
the construction of a connector which may be used with 
1852 (RG 114/U) cable. Measurements are also reported 
for 125Q2 attenuators and a 125Q pulse inverting transformer. 
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Theory of boil-off calorimetry, R. B. Jacobs, Rev. Sci. 
35, No. 7, 828-832 (July 1964). 

The usual formulation used to compute heat flow in boil-off 
calorimetry ignores the sensible heating of the vapor in the 
calorimeter. A formulation which includes this heating is 
derived; a single temperature measurement (of the vapor 
leaving the calorimeter) permits this effect to be computed. 
As the required temperature measurement has not been given 
for reported experiments, the vapor heating effect cannot be 
computed directly for those cases. Therefore, an expression 
for the computation of this temperature is derived theo- 
retically and applied to reported experiments. It is shown 
that, even in optimally designed calorimeters, the sensible 
heating of the vapor can be appreciable, and that the use of 
the usual formulation can introduce significant error in the 
computed heat flow. Finally, it is indicated how this error 
can explain, at least partially, the Heat Paradox which has 
been observed in insulation and Dewar evaluation tests. 


Instr. 


“Thumbtack”’ accelerometer for the 1.5-150 ke range, P. 


Hertelendy and P. Reynard, 
1305-1306 (Oct. 1964). — 

A special purpose mechanical vibration detector using a 
piezoelectric sensing element is described. Also, one appli- 
cation, a method of calibration, and some limitations of the 
device are discussed. Its chief merits are (1) low impedance, 
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(2) a nearly point contact, and (3) low transverse-to-longi- 
tudinal sensitivity. 


Microwave spectroscopy, |). | 
Phy. Chem. 15, 225-250 (1984). 
This article is a critical review of research in microwave 
spectroscopy of gasses during the past several years. The 
major topics covered are accurate structure determinations, 
structures of new or unusual molecules, internal rotation, 
vibrational information, high temperature spectroscopy, and 
instrumentation. A bibliography of the most important 
publications is included. 


Lide, Jr., 


Annual Rev. of 


The application of Archimedes’ principle to determination of 
the apparent volume of leather, J. R. Kanagy, J. Am. Leather 
Chem. Assoc. LEX, No. 11, 636-649 (Nov. 1964 

Studies were made of the volume changes occurring when 
acetone-dehydrated hide, chrome-retanned leather, vegetable- 
tanned strap leather, and chrome-retanned packing leathers 
were immersed in water, hydrocarbon solvent, 6M lithium 
bromide, and SM urea. None of these samples showed any 
volume changes in hydrocarbon solvent. All of the samples 
swelled considerably in water. All except vegetable-tanned 
strap leather showed shrinkage in 6M lithium bromide, and 
all showed some swelling in SM urea. 

Determinations of volume were made on specimens of hide 
and leather by immersing them in water after treatment with 
paraffin oil or dibutyl phthalate to fill the pores. These 
values were compared to values obtained for volumes by 
immersion in hydrocarbon solvent and by measurement of 
the dimensions of the specimens. The buoyancy methods 
were found to be quick, simple, and accurate procedures for 
determining volumes. Determinations in hydrocarbon sol- 
vent were found to average consistently approximately 0.8% 
higher than those obtained by filling the pores and immersing 
in water. 


Spectral tristimulus values for the CIE (u, v, w 
spacing system, I. Nimeroff, J. Opt. Soc. Am. 
54, No. 11, 1365-1367 (Nov. 1964 

The CIE uniform chromaticity space, u, v, w, developed by 
MacAdam and adopted in 1959, is known only by its chro- 
maticity transformations. The spacing should be known by 
its spectral tristimulus values as well if it is to be of greatest 
use in colorimetry. The set of trisimulus values, %, 2, wa, 
is found to have two bimodal functions, % and W. AS the 
CIE recommends, since 1963, the use of a 10° observer system 
as well as the 1931, 2° system, spectral tristimulus values 
are given for both systems. Also reported are the chroma- 
ticity coordinates of several standard sources, the locus of 
Planckian radiators, and the spectrum locus. 


uniform 
Le tter to Editor 


Change in the heat capacity of boron trioxide during the glass 
transformation, S. S. Chang and A. B. Bestul, J. Phys. Chem. 
68, 3082 (1964). 

A previous compilation showed that for B,O; the change in heat 
capacity during the glass transformation, Ac,, is only about 
one fifth of the ‘universal’? value which applies to most 
glass-forming substances. A computational error has been 
found in the value for B,O; in the above-mentioned compila- 
tion. The correct value of Ac, for BO; is 2.65 cal. deg. 
per “bead’’, which agrees very. well with Wunderlich’s 
‘“‘universal”’ Ac, value of 2.7+0.5 eal. deg.—'. Thus BO, in 
fact behaves like other glass-forming substances in this respect, 
rather than being exceptional. 


1 


Electrochemical calorimetry III. 
cadmium batteries, W. H 
Electrochem. Technology 2, 
1964). 

This paper reports some interesting sidelights obtained while 
making calorimetric measurements during various charges 
discharge cycles of nickel-cadmium cells proposed as satellite 
power sources by the National Aeronautics and Space 
Administration (NASA). Typical curves of heat changes as 
a function of time during operation are shown and discussed. 
The heat of reaction of oxygen with the constituents of the 
cell was measured and close agreement was obtained with the 
literature value for the heat of formation of cadmium oxide. 


Thermal effects of nickel- 
Metzger, Jr., and J. M. Sherfey, 
No. 9-10, 285-289 (Sept.—Oct. 
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Electrochemical calorimetry IV. ‘Measurement of the Peltier 
heat and voltage attending the passage of current across a 
liquid junction, A. Brenner and 8S. L. Gray, J. Electrochem. 
Soc. 111, No. 11, 1225-1230 (Nov. 1964). 

The experiments involved the measurement of the tempera- 
ture change resulting from the passage of current across the 
interface of two electrolytes. Measurements were made at 
various distances from the junction so as to vield a tempera- 
ture profile of the effect. The total heat effect was calculated 
from the integrated area under the curve of temperature 
profile. The coefficient obtained by dividing the total heat 
effect in joules by the current in coulombs was reported as a 
voltage equivalent. 

The following values were obtained for three systems of 
electrolytes and are probably correct to within about 7 mv: 


0.04 M hydrochloric acid and 0.4 M hydrochloric acid, 30 mv. 
0.90 M hydrochloric acid and 3.0 M potassium chloride, 80 


1.15 M hydrochloric acid and 1.21 M lithium chloride, 80 mv. 


\ hypothesis was advanced to relate the voltage equivalent 
of the heat effects to the diffusion liquid junction potentials 
and to the Peltier effect. 


The National Bureau of Standards tritiated water standards, 
W. B. Mann, J. 1 ppl. Rad. and Isotopes 15, 628, Pergamon 
Press Ltd. (1964) 

This letter notes the very agreement obtained by 
Euratom and the National Bureau of Standards in their 
calibration, by gas counting, of the National Bureau of 
Standards tritiated water standard. 


close 


Microwave and High-Frequency Calibration Services of the 
National Bureau of Standards, W. EF. Larson, EEE Trans. 
Microwave Theory Tech. MTT-12, No. 4, 480-482 (July 
1964) 

This is the second of a series on Microwave and High Fre- 
quency Calibration Services available at the National Bureau 
of Standards Boulder Laboratories, Boulder, Colorado. The 
quantities, attenuation, and field strength are described. 


Diffusion in a polymer with lamellar morphology, polyethyl- 
ene, R. K. Eby, J. Appl. Phys. 35, No. 9, 2720-2724 (Sept. 
1964). 

Data for the diffusion of ethane in polyethylene are presented 
together with the results of optical microscopy and of large- 
and small-angle x-ray diffraction for the same samples. 
The orientation of the lamella boundaries with respect to the 
concentration gradient of the diffusing molecules is important 
in determining the diffusion rate. In unoriented, as-crystal- 
lized polymers, the diffusion constant varies inversely with 
the lamella thickness. The results support the concept of 
lamella boundary diffusion. Lamella_ perfection is also 
important. In both oriented and unoriented polymers, the 
solubility constant apparently varies inversely with the 
lamella thickness, suggesting absorption in the lamella 
boundaries. 


Determination of the propagation constant for ultrasonic 
waves in melting and molten polyethylene, R. K. Eby, /. 
Acoust. Soc. Am. 36, No. 8, 1485-1487 (Aug. 1964). 

An extension of MeSkimin’s technique to polymer melts and 
soft amorphous polymers is reported. For 12 Me longitudinal 
waves, the modulus and internal friction are presented for 
the melting, supercooling, and annealing or recrystallization 
of polyethylene. Maxima of internal friction are observed 
in the solid at 52°, +60°, and between 100° and 130°C. 
The last, which is associated with the presence of crystals, 
decreases sharply near the melting temperature and is not 
present in the supercooled melt. In the melt, the internal 
friction exhibits a low, nearly constant value between 140 
and 290°C. In this range the ratio of the thermal coefficient 
of velocity change to the thermal coefficient of volume ex- 
pansion is 2.6. Above 300°C, there is a decrease of 
modulus and increase of internal friction believed to be 
associated with thermal degradation of the polymer. 


Lasers for length measurement, A. G. McNish, Science 146, 
No. 3541, 177-182 (Oct. 1964). 

CW lasers will find extensive use in length measurement 
because of the high coherence and intensity of the light they 
produce. ‘ 


Nonspattering solder flux, W. R. 
35, No. 10, 1369-1370 (Oct. 1964 
This article gives uses, methods of preparation and advantages 
of a polyethylene glycol soldering flux. 


Bjorklund, Rev. Sci. Instr. 


Simple are devices for spectral excitation in controlled 
atmospheres, M. Margoshes and B. F. Scribner, A ppl. Spectry. 
18, No. 5, 154-155 (1964). 

A gas-sheathed are is described for the excitation of spectra 
of solid materials in an atmosphere of any desired composition 
without the disadvantages associated with the chambers 
usually employed for this purpose. The atmosphere is con- 
trolled by enclosing the electrodes and the exeitation zone in 
a tubular flow of gas which effectively excludes air. The 
sample is held in an anode cap electrode whick can be changed 
rapidly between exposures. The are may be ignited within 
a few seconds after the flow of gas is started. The tubular 
flow of gas also stabilizes the are to a marked degree. 


Making precision measurements of zener diode voltages, 
W. G. Eicke, Jr., JEEE Trans. Commun. Elec. 83, No. 74, 
433-438 (Sept. 1964). 

The problems of making precise measurements of zener 
diode voltages are reviewed and discussed. Power supplies 
for the zener diode, measurement of the diode current, 
grounding of the circuit, shielding, mounting of the diode, 
and temperature control of the diode are discussed as they 
affect the measurement of the diode voltage. Equations are 
given that relate variations in diode voltage to variations in 
supply voltage and other parameters that affect the diode 
voltage. Three general classifications of methods for measur- 
ing diode voltages are defined. Specific techniques using equip- 
ment available in standards laboratories are described and 
equations for calculating the zener voltage and variation of 
zener voltage with variations in circuit parameters are given. 


Ethylene-propylene copolymers: crystallinity, infrared and 
creep studies, F. J. Linnig, EF. Parks, and L. A. Wood, J. 
Appl. Polymer Sci. 8, 2645-2651 (1964 

Density, infrared, and creep studies have been made of ethyl- 
ene-propylene copolymers containing up to about 50 mole 
“¢ propylene and their homopolymers. The degree of crystal- 
linity (volume fraction of polyethylene crystals), as measured 
by density and X-ray studies, decreases with increased propyl- 
ene content. This change is also reflected in the intensity 
of the 13.7 w band in the infrared spectra of the polymers. 
The measured amount of unsaturation (0.1-0.2° .) in the 
copolymers is too small to account for the presence of several 
infrared bands which are often attributed to unsaturation. 
An unsaturation of 0.90 was obtained for a terpolymer where 
the third component was unsaturated. Indentation-creep 
studies showed the creep to be nearly proportional to the 
logarithm of time with a slope of 10°. creep per decade for 
0.1 to 10 minutes. Compliance-temperature graphs were 
found to be quite similar to those for the styrene-butadiene 
rubber SBR 1500. 


The accurate measurement of voltage ratios of inductive 
voltage dividers, Acta /MEKO, pp. 317-331 (1964). Several 
methods have been developed at the National Bureau of 
Standards for the very accurate measurement of audio fre- 
quency voltage ratios associated with inductive voltage di- 
viders. Uncertainties in voltage ratio of less than 2 10-7 
have been obtained in measurements with transformer-type 
capacitance bridges. The attainment of such accuracy has 
provided a firm base for further investigations on inductive 
voltage dividers. An inherent major error arises from the 
interaction between distributed shunt impedances and leakage 
impedances in the windings. In specific designs of inductive 
voltage dividers the relative magnitudes of such errors are 
calculable and are characteristic of the design. A divider 
may be completely characterized by the designation of the 
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characteristic relative errors combined with an absolute | 
measurement. A method involving complementary measure- 
ments has been developed to yield an absolute measurement 
of the in-phase and quadrature components of the voltage 
ratio characteristic. Single-decade inductive voltage dividers 
have been constructed with ratio magnitude characteristics 
and residual errors less than 5X 10-%. This recent work has 
verified the accuracy of measurements made by other methods. 


Sampling and statistical design, W. J. Youden (Proc. 
Environmental Measurements, Cincinnati, Ohio, Sept. 4, 
Public Health Service No. 999-AP-15, pp. 35-39 (July 
A sampling scheme specifies the type of sample as well 
number and location of the samples. The placement 
samples depends on the region of interest. The 
sample depends on the accuracy or cost of the 
work. Composite samples may be used to conserve analytical 
effort. Alternatively, if the analytical errors are large, each 
sample should be large enough to permit multiple determina- 
tions. In any event the sampling scheme should be well 
adapted to statistical evaluation. For example, in order to 
detect trends in time or distance it is well to arrange for pairs 
of samples as follows 


Symp. 
1964), 
1964). 
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The paired dots represent samples taken at nearly the same 
time, or place, and the separation of the pairs corresponds to 
intervals in time or location. The ability to demonstrate 
heterogeneity in the region (temporal, spatial, or both) de- 
pends upon the agreement of adjacent samples. The pairs 
are not necessarily arranged in a straight line. Other arrange- 
ments will be discussed in detail. 


The needs of the statistical engineering programs of the 
National Bureau of Standards, C. Kisenhart, Proc. 9th Conf. 
Design of Experiments in Army Research Development and 
Testing, Directorate of Research and Development, U.S. Army 
Missile Command, Redstone Arsenal, Ala., Oct. 23 1963, 
ARO-D Report 64-2 (1964). 

The Statistical 


25, 


Engineering Laboratory of the Applied 
Mathematics Division is the largest of the three statistical 
engineering groups of the National Bureau of Standards. 
The staffing needs and recruitment difficulties of the Statisti- 
cal Engineering Laboratory are, it is believed, symptomatic 
of the problems faced by many statistical advisory and con- 
sulting groups in Government laboratories. It is pointed out 
that these recruitment difficulties stem in large part from 
the strong emphasis on theory and comparative neglect of 
applications in many undergraduate and_ postgraduate 
courses in mathematical statistics in the U.S.A., and from 
the absence in American colleges and universities of under- 
graduate programs leading to a bachelor’s degree in Statistics. 
The present recruitment problems of statistical engineering 
programs in Government laboratories are much as would be 
those of the corresponding physics programs if the only 
available physicists were MA’s and PhD’s in Theoretical 
Physies. 


Studies on the tungsten-rhenium thermocouple to 2000 °C, 
D. B. Thomas, Proc. 18th Annual ISA Conf. and Exhibit, 
Sept. 9-12, 1963, Chicago, Ill., Pre print 57.3.63, pp. 57.3.63 
1—57.3.63-9 (19683). 

Various lots of tungsten and rhenium wire were obtained 
from leading American manufacturers. Eleven tungsten- 
rhenium thermocouples were fabricated from these wire lots 
and tested at approximately 100 degree intervals up to 
2000 °C. Two polynomial equations were derived to rep- 
resent the average emf of the test thermocouples. Tables 
of temperatures versus emf are presented in 5 degree C 
intervals from 0 to 2000 °C and in 10 °F intervals from 32 to 
3640 °F. Spectrochemical analyses” of the tungsten and 
rhenium elements used in this study™are listed. A graphic 





comparison is made between the NBS emf values and values 
by other investigators. 


Modulus of rupture of glass in relation to fracture pattern, 
M. J. Kerper and T. G. Seuderi, Am. Ceramic Soc. Bull. 43, 


No. 9, 622-625 (Sept. 7, 1964). 
Glass which is broken in flexure shows a fracture face which 
normally consists in part of a smooth area which includes the 
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fracture origin. The area of this smooth portion or “mirror’ 
bears a relationship to the magnitude of the stress causing 
rupture. This relationship has been studied for eight different 
glasses tested under a variety of conditions. It was observed 
that the conditions of the glass surface, the test temperature, 
difference in exposure time to test temperature, and the 
rate of loading, all had an effect on the strength of glass but 
did not affect the relationship between modulus of rupture 
and the mirror size. The composition did, however, affect this 
relationship. This relationship is generally expressed as a 
power function where the modulus of rupture times the radius 
of the mirror to the one-half power is equal to a constant. 
This expression is not, however, valid for the full range of the 
data. 


Effect of particle shape and size distribution in a reinforced 
polymer, R. L. Bowen, J. Am. Dental Assoc. 69, 481-4965 
(Oct. 1964). 

Vinylsilane-treated silica-glass powders, bonded together with 
a self-curing monomer system, gave composites that could be 
used to restore anterior teeth. Spherical reinforcement 
particles required six to eight percent less resin (by volume) 
as a binder than did irregular-shaped particles. Intermittent 
or gap grading in the size distribution (three parts of fine 
combined with seven parts of coarse particles, medium sizes 
being omitted) eliminated an additional six to eight volume 
percent of the resin binder. Such a bimodal size distribution 
of spherical particles required only 13 percent resin (24 
volume percent); this composite had a relatively low coeffi- 
cient of thermal expansion and a high modulus of elasticity. 
The feasibility of using such composites (with reduced resin 
content) for dental restorations was supported, but the use of 
smaller particles was indicated. 

The shear strength of human dentin was measured, using a 
punch test on transverse sections of tooth specimens. The 
average shear strength was about 14,000 pounds per square 
inch. 


Free-field correction for condenser microphones, W. Koidan 
and D. 8. Siegel, J. Acoust. Soc. Am. 36, No. 11, 2233 
(Nov. 1964). 

Differences in the free-field corrections of one-inch condenser 
microphones, caused by their different acoustic driving- 
point impedances, were found to be of sufficient magnitude to 
warrant a cautious approach to the determination of the 
free-field response level from measurements of the pressure 
response level and the addition of a “‘standard’”’ free-field 
correction. 
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A magnetic fixed point for thermometry below]l °K, J. C. 
Kisenstein, R. P. Hudson, and B. W. Mangum,*A ppl. Phys. 
Letters 5, No. 11, 231 (Dec. 1, 1964). 

It is suggested that a peak in the magnetic susceptibility of 
CeClh, at 0.345 °K would be useful as a thermometric fixed 
point. 


Some activities of the cooperative dental research between the 
Federal Government and the American Dental Association 
at the National Bureau of Standards, G. C. Paffenbarger, 
Ann. Dentistry XIII, No. 4, 111-126 (Dec. 1964). 

Some aspects of the cooperative dental research between the 
Federal Government and the American Dental Association 
conducted at the National Bureau of Standards inelude 
(1) the necessary research and testing to formulate specifica- 
tions for dental materials, (2) the establishment of a List of 
Certified Dental Materials by the Association, (3) the develop- 
ment of diagnostic instruments, (4) the production of motion 
pictures depicting how the research findings may be applied 
in practice, (5) the determination of physical and chemical 
properties of dental materials and the hard tooth tissues, 
(6) the development of new, and improvements in old, 
materials, and (7) establishment of the validity of some labora- 
tory findings. Examples in each category are presented. 


Studies of the thicknesses of adsorbed} glass finished by 
ellipsometry, D. J. Tutas, R. R. Stromberg, and E. Passaglia, 
SPE Trans. 4, No. 4, 1-7 (Oct. 1964). 





The thickness of the adsorbed layer of two glass finishes, 
vinyl tris(2-methoxyethoxy)silane and gamma aminopropyl- 
triethoxysilane, was measured by ellipsometry. The vinyl 
silane was adsorbed on glass from solutions in water and in 
methyl ethyl ketone. The adsorption from water resulted in 
a film initially at least 50 A thick, which increased in thickness 
with time. Adsorption from methyl ethyl ketone solution 
resulted film of about 10 A. The 
sorbed on and on chrome from 


amino silane was ad- 
solutions in water. 
A in thick- 
remained constant over a period of about 5 hours. 
Adsorption on chrome resulted in a film about 10 A in thick- 
From the results it appears that a polymer, rather than 
the small molecules of glass finish, are adsorbed from the 
water solution on glass and that small unpolymerized ‘“‘finish 
molecules are adsorbed from the methyl ethyl ketone solution. 


glass 


Adsorption on glass resulted in a film at least 60 
ness, which 


ness. 


Influence of “or energy sources on bitumen oxidation, 
Kk. G. Martin, Appl. Chem. 14, No. 11, 514 (1964). 

Using the cae ‘change in infrared absorption of the carbonyl 
functional groups of asphalt with increase of radiant energy 
flux, the oxidation rates of four air-blown asphalts of coating 
grade have been obtained for the three sources of radiation, 
carbon are, xenon are, and sunlight. Differences in the 
pattern of oxidation for the much used carbon are compared 
to both xenon-are and solar irradiation were found, whereas 
the latter two sources induced nearly identical behaviour in 
the asphalts. Continuous measurement of the intensity of 
the solar radiation by pyroheliometry has indicated that the 
asphalt oxidation rates are directly proportional to intensity, 
sensitive to very low intensities and independent of diurnal 
changes. Potassium ferrioxalate dosimetry has shown that 
the intensity of the xenon-arc irradiation as used in a weather- 
simulated that of summer sunshine extremely well, 
whereas the carbon-are irradiation as used in the ASTM 
method of evaluating weatherometer durability of asphalts 
had a fivefold greater intensit) This increased intensity 
changed the pattern of rate of asphalt oxidation rather than 
simply accelerating the steady-state oxidation rates. 


ometer 


Which measure of precision? The evaluation of the precision 
of ‘Cie "r= involving linear calibration curves, 
F. Linnig and J. Mandel, Anal. Chem. 36, No. 13, 25 A—-32. 
ety 1964 

Three sources of error can affect 


the precision of an analytical 
procedure involving calibration 


lines: (1) variability among 
replicates, (2) seatter about the calibration line, and (3 
uncertainty of the calibration line In many cases, the first 
of these is small in comparison with the other two; yet, it is 
often used as the measure of the precision of the 


analytical 
method. A detailed analysis is given of the three 


sources of 


error, and methods of calculation are presented for the segre- 
gation of these errors on the basis of a calibration experiment. 


Other NBS Publications 


Journal of Research 69A 
Apr. 1965), 70 cents. 
tadiolysis of NUN4O. R. Gorden, i and P. Ausloos. 
Oscillator strengths for lines of Ni 1 . H. Corliss. 
More — ‘lengths from thorium vst W. F. 
R. W. Stanley. (See above 
“Fall equilibria in the system vanadium 
oxide. J. L. Waring and R. SS. Roth. 
Thermodynamics of oe? ternary system: 
sium chloride at 25 
. E. Bower and R. Robinson. 
Heat of formation of ee fluoride by 
aged elements. E.S i 
Relative enthal; Vv of polytet rafluoro¢ thyle ne 
C. T. B. Douglas and A. W. Harman. 
Anionic polymerization of isoprene at low 
polyisoprenyllithium. L. J 
One particle transitions 
chanies. A. R. Ruffa. 
Disaccommodation of magnetic spectra of two 
zine ferrites. A. L. Rasmussen. 


Phys. and Chem.), No. 2 (Mar.- 
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—— vapor pressure measurements. 


direct 
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Fetters. 
and correlation in 


quantum me- 


manganese 


Splitting of a set of equivalent sites in centrosymmetric 
space groups into subsets under homogeneous stress. J. B 
Wachtman, Jr., and H. 8. Peiser. (See above abstracts.) 


Radio Sci. 

$1.00. 

Klectromagnetic wave penetration of reentry plasma sheaths. 
M. P. Bachynski. 

On the use of refractive index diagrams for 
anisotropic regions. L. B. Felsen. 

Surface waves along a perfectly conducting plane 
with semi-infinite magneto-plasma. 8S. Adachi 
Mushiake. 

On the point of emergence of a microwave 
linearly graded plasma. <A. L. Cullen. 

Self and mutual admittances of waveguides 
plasma layers. J. Galejs. above 

kffect of electron collisions on the 
theory. K.G. Budden. 

Momentum transfer collisions in 
trons. M. H. Mentzoni. 

I:xperimental studies of perturbations in ionospheric plasma. 

H. Heisler. 

Electromagnetic 

axial magnetic 
H. Walter. 

Radar sections of 
quency. L. Peters, Jr 

Klectroacoustic waves excited by a space 
atmosphere and its effect on radar 
(See above abstracts. 

Discussion on basic equations with 
pressive plasmas. Kk.-M. Chen 

Scattering of electromagnetic and electroacoustic 
a evlindrical object in a compressible plasma. J. R. Wait. 

Multiphase periodic very-low-frequency emissions. N. Brice. 

Radio studies of the high-latitude ionosphere during the solar 
eclipse of 20 July 1963. R. D. Hunsucker. 

VLF and LF fields propagating near and into a rough sea. 
R. M. Lerner and Max 

Insulated and loaded loop antenna immersed in a conducting 
medium. R.H. Williams. (See above abstracts.) 

‘apacitance of biconical antennas in magneto-ionic 
elliptic cone capacitance. V. P. Pyati and H. Weil. 

‘alculations of the bistatic scattering section of 
sphere with an impedance boundary condition. J. 
Wait and C. M. Jackson. 

‘oncept of differential reflectivity as applied to the reflection 
of beam-limited radiation by a convex body. A. Erteza, 
J. A. Doran, and D. H. Lenhert. 


J. Res. NBS/USNC-URSI 69D, No. 2 (Feb. 1965), 
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Radio Sci. J. Res. NBS/USNC 
$1.00. 

Propagation in nonuniform gyrotropic 
and L. B. Felsen. 
Geometrical optics for 
L. Peters, Jr 
Attenuation of 

S. I. Akasofu. 

Self distortion of radio signals in the 
Atmospheric gravity waves: 
C. O. Hines. 
Electromagnetic 
R. KE. Burgess. 
A generalized hydromagnetic wave in an inhomogeneous, 
evlindrical plasma. C. K. McLane and T. Tsukishima. 
Angular dependence of the refractive index in the ionosphere. 

G. A. Deschamps. 
Electrodynamics of moving anisotropic media: the 
theory. CO, 2. cai. 
Study of the phenomenon of whistler echoes. 
W. C. Johnson, and J. F. Walkup. 
Multiple-frequency investigations of radio wave 
during the dawn-breakup phase of auroras. R. 
rathy, and F. T. Berkey. 
Sferiec excitation of a two-laye! 
Kraichman. 


URSI 69D, No. 3 (Mar. 1965), 


media. S. H. Gross 


gyrotropic bodies. W. ;  s Lee, 
, and C. H. Walter. 
hydromagnetic waves in the ionosphere. 
D region. L.R 
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Laaspere, 


absorption 
Parthasa- 


conducting medium. M. B. 


160 





Currents, charges, and near fields of cylindrical antennas. 
R. W. P. King and T. T. Wu. 

A note on the radiation conductance of an axial slot on a 
evlinder. C. M. Knop and C. T. Swift. 

Analytical formulas for radio paths in spherically stratified 
ionospheres. Ek. Woyk (Chvojkova 


Radio Sci. J. Res. NBS/USNC-URSI 69D, No. 4 (Apr. 1965), 
$1.00. 


VW aves in Plasma Pape rs 


Dispersion of waves in a cold magneto-plasma from hydro- 
magnetic to whistler frequencies. H. G. Booker and R. B. 
Dyce. 

Study on the guiding mechanism of whistler radio waves. 
S. Adachi. 

Electromagnetic waves along an infinitely long and thin con- 
ducting wire in a magneto-ionie medium. Y. Mushiake. 
(See above abstracts. 

Use of the phase-integral method to determine the reflection 
properties of a stratified ionosphere. C. Altman. 

E-mode propagation in a plane-stratified plasma. 
and J. Shmoys 

fadiation from an infinite axial slot on a circular cylinder 
clad with magnetoplasma. P. de Marchin and G. Tyras. 

Index of refraction surfaces for plasma waves. T. Yeh and 
M. HH. Cohen. 

Impedance of a short dipole in a compressible plasma. K. 
G. Balmain 

Waves circulating around a rigid cylindrical obstacle in a 
compressible plasma. J. R. Wait. 

Wave propagation in a two component warm plasma. 
Seshadri. 


P. Hirseh 


Harmonie currents excited by an electromagnetic wave in a 
plasma. L. Wetzel and T. W. Tang. 

Excitation of acoustic waves in plasmas. W. 
(See above abstracts.) 

Wave interaction in 
Narasinga Rao. 


A. Saxton. 


oxygen magneto plasmas. K. 


Regular Contributions 


Test of the constancy of the velocity of electromagnetic radia- 
tion in high vacuum. P. Beckmann and P. Mandies. 
(See above abstracts). 

Signal degeneration in laser beams propagated through a 
turbulent atmosphere. P. Beekmann. 
stracts). 

High frequency backscatter from the earth measured at LOOO 
km altitude. R. C. Chia, A. K. Fung, and R. K. Moore. 

Observation of NPG VLF transmissions at Tracy, California 
during path equinox. G. B. Carpenter and A. L. Whitson. 

Small magnetic toroid antenna imbedded in a highly con- 
ducting half space. G. R. Swain. ; 


(See above ab- 


Radio Sci. J. Res. NBS/USNC-URSI 69D, No. 5 (May 1965), 

$1.00. 

Preface to fourth of “Waves in Plasma”’ papers. J. R. Wait 

On electromagnetic radiation from a magnetic dipole with 
arbitrary orientation embedded in a lossless magneto-ionic 
medium. H. Motz. 

A systematic study of the radiation patterns of a dipole in 
magnetoplasma based on a classification of the associated 
dispersion surfaces. R. Mittra and G. L. Duff. 

Propagation of vertically polarized electromagnetic waves in 
a horizontally stratified magnetoplasma. R. Burman and 
R. N. Gould. 

A note concerning the reflection of waves in inhomogeneous 
lavers with asymmetric profiles. R. Burman. 

On the terrestrial propagation of ELF and VLF waves in the 
presence of a radial magnetic field. J. Galejs. 

Electron density profiles in cylindrical plasmas from micro- 
wave refraction data. B. A. Aniéin. 

Study of the waves supported by a warm plasma slab. P. R. 
Caron. 


Transmission and reflection of electromagnetic waves by a hot 
plasma. KE. C. Taylor. 

Radiation from electrons in a 
Liemohn. 

fadiation from a uniformly moving charge in an anisotropic, 
two component plasma. 8S. R. Seshadri and H. 8. Tuan. 


magnetoplasma. H. B. 


Microwave spectral tables line strengths of asymmetric rotors, 
P. F. Wacker and M. R. Pratto, NBS Mono. 70, Vol. II 
(Dee. 15, 1964), $3.00. 

On the statistical theory of electromagnetic waves in a 
fluctuating medium (II). Mathematical basis of the 
analogies to quantum field theory, K. Furutsu, NBS Mono. 
79 (Dee. 7, 1964), 35 cents. 
complete mode sum for LF, VLF, ELF 
wave fields, J. R. Johler and L. A. Berry, 
(Oct. 1, 1964), 20 cents. 

A study of lunar surface radio communication, L. E. 
NBS Mono. 85 (Sept. 14, 1964), 70 cents. 

Standard Reference Materials: Preparation of NBS copper- 
base spectrochemical standards, R. KE. Michaelis, L. L. 
Wyman, and R. Flitsch, NBS Mise. Publ. 260-2 (Oct. 15, 
1964), 35 cents. 

Quarterly radio noise data June, July, August, 1963, W. Q. 
Crichlow, R. T. Disney, and M. A. Jenkins, NBS Tech. 
Note 18-19 (Aug. 21, 1964), 50 cents. 

Quarterly radio noise data September, October, November 
1963, W. Q. Crichlow, R. T. Disney, and M. A. Jenkins, 
NBS Tech. Note 18-20 (Oct. 23, 1964), 50 cents. 

Calculation of sunrise and sunset times at ionospheric heights 
along a great circle path, A. H. Brady and D. D. Crombie, 
NBS Tech. Note 209 (Nov. 8, 1964), 20 cents. 

Demagnetizing factors for oblate spheroids used in ferri- 
magnetic resonance Measurements, L. B. Schmidt, W. E. 
Case, and R. 1D. Harrington, NBS Tech. Note 221 (Sept. 4, 
1964), 20 cents. 

Concerning the theory of radiation from a slotted conducting 
Plane in a plasma environment, J. R. Wait, NBS Tech. 
Note 223 (Sept. 28, 1964) 

An atlas of VLF emission spectra observed with the ‘Hiss 
recorder,’ J. A. Koch and V. C. Edens, NBS Tech. Note 
226 (Nov. 13, 1964), 40 cents. 

VHF ionospheric scatter system loss measurements European- 
Mediterranean area, V. H. Goerke and O. D. Remmler, 
NBS Tech. Note 230 (Dee. 25, 1964), 40 cents. 

A program for plotting circles of constant overpressure around 
targeted points, M. L. Joel and D. D. Lottridge, NBS 
Tech. Note 249 (Oct. 28, 1964 $0 cents. 
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(Oct. 2, 1964), 15 cents. 

On the pictorial structure of Chinese characters, B. K. Rankin 
III, W. A. Sillars, and R. W. Hsu, NBS Tech. Note 254 
Jan. 4, 1965), 30 cents. 

Characteristics of the earth-ionosphere waveguide for VLF 
radio waves, J. R. Wait and K. P. Spies, NBS Tech. Note 
300 (Dee. 30, 1964), 50 cents. 

Interpretation of potentiometric titrations of weak acids in 
methanol-water solvents, K. C. Ong, R. A. Robinson, 
and R. G. Bates, Anal. Chem. 36, 1971-1972 (Sept. 1964 

Fluctuations in ionospheric absorption events at conjugate 
stations, J. K. Hargreaves and H. J. A. Chivers, Nature 
203, No. 4948, 963-964 (Aug. 29, 1964 

RF and microwave power measurements, G. F. Engen, Proce. 
URSTI XIVth General Assembly, Tokyo, Japan, Sept. 1963, 
XIII, 77-80 (1964). 

LF—-VLF frequeney and time services of the National Bureau 
of Standards, ID. H. Andrews (Proce. 19th Annual ISA 
Conf. and Exhibit, New York, N.Y., Oct. 12-15, 1964 
ISA Preprint 21. 3—-2—64 (1964). 

Fresnel zone diffraction effeets at 50 Ge/see. determined from 
measured aperture field data, R. C. Baird, Program and 
Digest 1963 IEEE Prof. Group Ant. and Prop. Symp. 


Bowman, NBS Tech. Note 250 





Space Telecommunications, Boulder, | 
Colo., pp. 171-173 (1963). 

Electronic energy bands in strontium titanate, A. 
and A. J. Leyendecker, Phys. Rev. 135, No. 5A, 

Aug. 31, 1964). 

On the implication of diurnal, seasonal and geographical 
variations in composition of high atmosphere from F- 
region measurements, J. W. Wright (Proc. NATO Ad- 
vanced Study Institute, Sheikampen, Norway, Apr. 17-26, 
1963), Book, Electron Density Distribution in the TIono- 
sphere and Exosphere, pp. 186-198 (North Holland Publ. 
Co., Amsterdam, The Netherlands, 1964). 

Ionospheric winds: motions into night and sporadic E cor- 
relations, N. W. Rosenberg, H. D. Edwards, and J. W. 
Wright (Proc. 4th Intern. Space Science Symp., Warsaw, 
June 4-10, 1963), Book, Space Research, ed., P. Muller, 
IV, 171-181 (North Holland Publ. Co., Amsterdam, The 
Netherlands, 1964). 

Oblique incidence pulse measurement at 
Dougherty, AGARDograph 74, 
Waves at Frequencies below 300 KC, p. 133 
gamon Press Inc., New York, N.Y., 1963). 

Nonresonant absorption and collision diameters in the foreign- 
gas broadening of symmetric top molecules, G. Birnbaum 
and A. A. Marvott, J. Chem. Phys. 41, No. 1, 154-157 
July 1, 1964). 

Electromagnetic cross sections for electron and nuclear re- 
search, H. W. Koch, Nucl. Instr. Methods 28, 199-204 
1964). 

Exclusion of parity unfavored transitions in forward scatter- 
ing collisions, U. Fano. Phys. Rev. 135, No. 3B, B863 
B864 (Aug. 10, 1964). 

Microscopic-astrophysics and the development of astro- 
physical-laboratories, R. N. Thomas, Bull. Centre Intern. 
d'Astrophysique del’Observatoire de Nice, No. 1, pp. 7-16 
Nice, France, 1964). 

Alouette ionic, magnetic-field and plasma studies, W. Calvert, 
H. Rishbeth, and T. E. VanZandt, IG Bull. No. 83, 16-19 
May 1964); Trans. AGU 45, No. 2, 398-401 (June 1964). 

Propagation of radio waves past a coast line with a gradual 
change of surface impedance, J. R. Wait and K. P. Spies, 
IEEE Trans. Ant. Prop. AP-12, No. 5, 570-575 (Sept. 
1964). 

On the propagation of 
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longwave-length terrestrial radio 
waves—two theoretical techniques, J. R. Johler and L. A. 
Berry, AGARDograph 74, Propagation of Radio Waves 
below 300 KC, pp. 101-129 (Pergamon Press Inc., New 
York, N.Y., 1963). 

Experiments on self-ionizing shock waves in a magnetic field, 
Y. Nakagawa and K. B. Earnshaw, Book, Comptes Rendus 
de la Vie Conference Internationale sur les Phenomenes 
D’Ionisation dans les Gaz, ed. P. Hubert and E. Cremieu- 
Alean, IV, 303-309 (S.E.R.M.A., Paris, France, 1963). 

Microwave propagation in an overdense bounded magneto- 
plasma, B. Wieder, Phys. Fluids 7, 964-972 (July 1964). 

Ionosonde studies of some chemical releases in the ionosphere, 
J. W. Wright (Proc. NATO Advanced Study Institute, 
Sheikampen, Norway, Apr. 17-26, 1963), Book, Electron 
Density Distribution in the Ionosphere and Exosphere, pp. 
116-137 (North Holland Publ. Co., Amsterdam, The 
Netherlands, 1964). 

Standard nuclear instrument modules, L. Costrell, U.S. 
Atomic Energy Commission TID Report 20893, 12 pages 
(July 1964 

Charge-transfer absorption spectra of NO in Kr and CH;OH 
solutions, E. E. Ferguson and H. P. Broida, J. Chem. 
Phys. 40, No. 12, 3715-3716 (June 15, 1964). 

Ionospheric storms, 8. Matsushita, Proc. URSI 
General Assembly. Tokyo, Japan, Sept. 1963, 
169-180 (1964). 

Symplectiec modulary groups, M. Newman and J. R. Smart, 
Acta Arithmetica IX, 83-89 (1964). 

Standards and measurement of attenuation, impedance and 
phase shift, R. W. Beatty, Proce. URSI XIVth General 
Assembly, Tovko, Japan, Sept. 1963, XIII, 80-92 (1964 

Vapour pressure, G. T. Armstrong, Eneyclopaedie Dictionary 
of Physics 7, 586-587 (Pergamon Press Ine., New York, 
N.Y., 1962). 
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Exterior corrosion of cast-iron pipe, M. Romanoff, J. Am. 
Water Works Assoc. 56, No. 9, 1129-1143 (Sept. 1964). 
The energy loss of electrons in solids, A. M. Rauth and J. A. 

Simpson, Radiation Res. 22, No. 4, 643-661 (Aug. 1964). 

Coherence theories of tropospheric radio propagation, J. 
Wait, IEEE Trans. Ant. Prop. AP-12, No. 5, 
(Sept. 1964). 

Solar activity, J. V. Lincoln, IQSY Instruction Manual. 
No. 2, Report of sudden ionospheric disturbances, Suppl., 
pp. 1-11 (IQSY Secretariat, London, England, 1963). 

Vibration-rotation bands of deuterium cyanide and hydrogen 
evanide, A. G. Maki, E. K. Plyler, and R. Thibault, J. 
Opt. Soc. Am. 54, No. 7, 869-876 (July 1964). 

IQSY Instruction Manual. No. 1. World days, A. H. 
Shapley, 32 pages, IQSY Secretariat, London, England 
(July 1963). 

Influence of the refractive index profile in VHF reflection 
from a tropospheric layer, J. R. Wait and C. M. Jackson, 
IEEE Trans. Ant. Prop. AP-12, No. 4, 512-513 (July 1964). 

Photometric observations of the airglow during the IQSY. 
F. E. Roach, IQSY Instruction Manual. No. 5: Air- 

pt. I, 1-24 (IQSY Secretariat, London, England, 
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1963). 

Correlation matrix for radio waves in the troposphere, 
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